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Palavras-chave Problemas variacionais fraccionais em tempo discrete, operadores frac- 

cionais discretes do tipo de Riemann-Liouville, soma por partes frac- 
cional, equacoes de Euler-Lagrange, condicao necessaria de optimali- 
dade do tipo de Legendre, escalas temporais, sistema de computacao 
algebrico Maxima. 

Resumo Introduzimos um calculo das variances fraccionais nas escalas tem- 

porais Z e (hZ)^. Estabelecemos a primeira e a segunda condicao 
necessaria de optimalidade. Sao dados alguns exemplos numericos que 
ilustram o uso quer da nova condicao de Euler-Lagrange quer da nova 
condi^iao do tipo de Legendre. Introduzimos tambem novas defini^ioes 
de derivada fraccional e de integral fraccional numa escala temporal 
com recurso a transformada inversa generalizada de Laplace. 
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Riemann-Liouville fractional-order operators, fractional formula for 
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Abstract We introduce a discrete-time fractional calculus of variations on the 

time scales Z and {hZ)a- First and second order necessary optimality 
conditions are established. Some numerical examples illustrating the 
use of the new Euler-Lagrange and Legendre type conditions are 
given. We also give new definitions of fractional derivatives and 
integrals on time scales via the inverse generalized Laplace transform. 
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Introduction 



The main goal of this thesis is to develop a more general fractional calculus on time 
scales. 

In the first year of my PhD Doctoral Programme I followed several one-semester courses 
in distinct fields of mathematics. One of the one-semester courses was called Research Lab 
where five different subjects were covered. One of those subjects was calculus of variations 
on time scales. Time scales theory looked to me a new and wonderful subject. As the 
doctoral programme PDMA Aveiro-Minho is backed up by two research units, not only 
Professors of the first semester but any researcher that belongs to any of these R&D units is 
invited to present, at the end of the first semester, research topics for PhD thesis by giving 
seminars to interested students. The author's supervisor was one of the researchers that 
was available to be advisor if any student show interest in the theme Fractional Calculus 
on Time Scales - the title of this thesis. 

The main idea was to connect in one theory two subjects that were, and still are, subject 
to strong research and development. Our contribute on this new area, until now, appears 
in part II of this thesis. On one hand we develop a discrete fractional calculus of variations 
for the time scale Z (Chapter [3]) and for the time scale {h'L)a (Chapter Hj). One the other 
hand we give new definitions for differintegral fractional calculus on a time scale using a 
Laplace transform approach (Chapter |5]). 

Nowadays fractional differentiation (differentiation of an arbitrary order) plays an im- 
portant role in various fields: physics (classic and quantum mechanics, thermodynamics, 
etc.), chemistry, biology, economics, engineering, signal and image processing, and control 
theory p[6ni[72l[73ll871IM[95l fTnnl[TM[TTn] . It is a subject as old as Calculus itself but much 
in progress |98]. The origin of fractional calculus goes back three centuries, when in 1695 
L'Hopital asked Leibniz what should be the meaning of a derivative of order 1/2. Leibniz's 
response: "An apparent paradox, from which one day useful consequences will be drawn". 
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After that episode, which most authors consider the born of fractional calculus, several fa- 
mous mathematicians contributed to the development of Fractional Calculus iGO l lHT tfTOS] : 
Abel, Fourier, Liouville, Griinwald, Letnikov, Caputo, Riemann, Riesz, just to mention a 
few names. In the last decades, considerable research has been done in fractional calculus. 
This is particularly true in the area of the calculus of variations, which is being subject to 
intense investigations during the last few years [IIl[25l|26l[92l[l02l[l03] . F. Riewe [TnKlfTnH] 
obtained a version of the Euler-Lagrange equations for problems of the calculus of varia- 
tions with fractional derivatives, that combines the conservative and non- conservative cases. 
In fractional calculus of variations our main interest is the study of necessary optimality 
conditions for fractional problems of the calculus of variations. The study of fractional 
problems of the calculus of variations and respective Euler-Lagrange equations is a fairly 
recent issue - see [SHHElElIIDlIIlllSSllinilSDllSSllSSllHI] and references therein - and include 
only the continuous case. It is well known that discrete analogues of differential equations 
can be very useful in applications [28 | [68 | [7T] and that fractional Euler-Lagrange differen- 
tial equations are extremely difficult to solve, being necessary to discretize them [6|I26]. 
Therefore, we consider pertinent to develop a fractional discrete-time theory of the calculus 
of variations in a different time scale than M. We dedicate two chapters to that: one for 
the time scale Z and another for the time scale {KL)a- Applications of fractional calculus 
of variations include fractional variational principles in mechanics and physics, quanti- 
zation, control theory, and description of conservative, nonconservative, and constrained 
systems [2 ^ 150 1 15T |I1U3] . Roughly speaking, the classical calculus of variations and optimal 
control are extended by substituting the usual derivatives of integer order by different kinds 
of fractional (non-integer) derivatives. It is important to note that the passage from the 
integer/classical differential calculus to the fractional one is not unique because we have 
at our disposal different notions of fractional derivatives. This is, as argued in |25 |I1U2| . an 
interesting and advantage feature of the area. Most part of investigations in the fractional 
variational calculus are based on the replacement of the classical derivatives by fractional 
derivatives in the sense of Riemann-Liouville, Caputo and Riesz [SlIinillSlESlEZlES]- In- 
dependently of the chosen fractional derivatives, one obtains, when the fractional order of 
differentiation tends to an integer order, the usual problems and results of the calculus of 
variations. 

A time scale is any nonempty closed subset of the real line. The theory of time scales 
is a fairly new area of research. It was introduced in Stefan Hilger's 1988 Ph.D. thesis |61] 
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and subsequent landmark papers [621163] . as a way to unify the seemingly disparate fields of 
discrete dynamical systems (i.e., difference equations) and continuous dynamical systems 
(i.e., differential equations). His dissertation referred to such unification as "Calculus on 
Measure Chains" [2^175]. Today it is better known as the time scale calculus. Since the 
nineties of XX century, the study of dynamic equations on time scales received a lot of 
attention (see, e.g., [2ll?llll3]). In 1997, the German mathematician Martin Bohner came 
across time scale calculus by chance, when he took up a position at the National University 
of Singapore. On the way from Singapore airport, a colleague, Ravi Agarwal, mentioned 
that time scale calculus might be the key to the problems that Bohner was investigating at 
that time. After that episode, time scale calculus became one of its main areas of research. 
To the reader who wants a gentle overview of time scales we advise to begin with that 
was written in a didactic way and at the same time points some possible applications (e.g. 
in biology). Here we are interested in the calculus of variations on time scales |85 [I112| . 
Our goal is to connect the theories of calculus of variations on time scales and fractional 
calculus. 

This thesis is divided in two major parts. The first part has two chapters in which we 
provide some preliminaries on fractional calculus and time scales calculus, respectively. The 
second part is splitted in three chapters where we present our original work. In Chapter [3] 
we introduce the fractional calculus of variations on the time scale Z. The main results 
of this chapter are a first-order necessary optimality condition (Euler-Lagrange equation) 
and a second-order necessary optimality condition (Legendre inequality). In Chapter H] we 
introduce a fractional factorial function that allow us to define left and right fractional 
derivatives and to develop further the two necessary optimality conditions of Chapter [3l 
In Chapter m the time scale is /iZ^. We believe that our results open some possible doors 
to research, also in the continuous case, when h ^ 0. Chapter O is devoted to a new 
approach to define fractional derivatives and fractional integrals in an arbitrary time scale, 
by using the Laplace transform of time scales as support. Finally, we write our conclusions 
in Chapter E] as well some future research directions. 
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Part I 
Synthesis 



Chapter 1 



Fractional Calculus 



"The fractional calculus is the calculus of the XXI-st century" 

K. Nishimoto (1989) 

The theory of discrete fractional calculus is in its infancy [20 | l2T | [86]. In contrast, the 
theory of continuous fractional calculus is much more developed |1U8| . The fractional 



discrete theory has its foundations in the pioneering work of Kuttner in 1957, where it 



appears the first definition of fractional order differences. In Sections 11.11 and 11.21 some 
definitions of fractional discrete and continuous operators are given, respectively. 

1.1 Discrete fractional calculus 

In 1957 [71] Kuttner defined, for any sequence of complex numbers, {an}, the s-th order 
difference as 

oo 



m=0 

where 



■^-^ f —s — 1 -\- m\ 

A On = ( jan+m, (1-1 

m=0 \ ^ / 

t\ t{t - 1) . . .{t - m+ 1) 



ml ml 



1.2) 



In [71] Kuttner also remarks that ^ ^ means when m is negative and also when 
t — m is a negative integer but t is not a negative integer. Clearly, (II. ip only makes sense 



when the series converges. 
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In 1974, Diaz and Osier |48] gave the following definition for a fractional difference of 
order 

" ^ " |/(x + z/-/c), 



; — n V / 



fc=0 



^kj T{u-k + l)k\ 
where u is any real or complex number. 

The above definition uses the usual well-known gamma function which is defined by 



:i.3) 



OO 



T{u) := / f'^'e^'dt, V e C\(-No) . (1.4) 
Jo 

The gamma function was first introduced by the Swiss mathematician Leonard Euler 
in his goal to generahze the factorial to non integer values. 

Throughout this thesis we use some of its most important properties: 

r(i) = 1, (1.5) 

r(z/ + l) = uT{u) for z/ e C\(-No), (1.6) 
r(x + 1) = x\, for X G No . (1.7) 

Remark 1. Using the properties of gamma function it is easily proved that formulas (II. 2p 
and (II. 3p coincide for u an integer. 

We begin by introducing some notation used throughout. Let a be an arbitrary real 
number and h = a + k for a certain A; G N with k > 2. Let T = {a, a + 1, . . . ,b}. According 
with mi, we define the factorial function 



=t(t-l)(t-2)...(t-n + l), neN 

and = 0. 

Extending the above definition from an integer n to an arbitrary real number a, we have 

. -pi±^, (1.8) 

r(t+l-a)' ^ ' 

where F is the Euler gamma function. In [86] Miller and Ross define a fractional sum of 
order u > via the solution of a linear difference equation. Namely, they present it as 
follows: 
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Definition 2. 

^ t—f 

^-"fit) = fZT - <s)r-'^f{s). (1.9) 

^ ' s=a 

Here / is defined for s = a mod (1) and A^*^/ is defined for t = {a + u) mod (1). 

Tliis was done in analogy witli tlie Riemann-Liouville fractional integral of order z/ > 
(cf. formula fll.lOp ) which can be obtained via the solution of a linear differential equation 
[86|[87]. Some basic properties of the sum in (11 ■9p were obtained in [86]. Although there 
are other definitions of fractional difference operators, throughout this thesis, we follow 
mostly, the spirit of Miller and Ross, Atici and Eloe PT|I86]. 



1.2 Continuous fractional calculus 

In the literature there are several definitions of fractional derivatives and fractional inte- 
grals (simply called differintegrals) like Riemann-Liouville, Caputo, Riesz and Hadamard. 
Throughout this thesis only Riemann-Liouville and Caputo definitions are used. 
By historical precedent over Caputo definition, we will start with the Riemann-Liouville 
definition. Let [a, h] be a finite interval and a G M"*". The left and the right Riemann- 
Liouville fractional integrals (RLFI) of order a of a function / are defined, respectively, 
by: 

ai:fix):=-^ [ {x-tr-'f{t)dt, x>a, (1.10) 

and ^ 

Jbfix) := FTT / (t- ^r~'fit)dt, x<b, (1.11) 

where T{a) is the Gamma function (II. 4p . 

Remark 3. If / is a continuous function al^f ~ x^t ~ f ■ 

Remark 4. The semigroup property of Riemann-Liouville fractional operators are given by 

{ai:al^f){x) = ai:^^f{x), X > a, « > 0, /3 > 0. (1.12) 

If / is a continuous function a^xf — ~ f ■ 

Let / be a function, a G and n = [a] + 1, where [a] means the integer part of a. 
The left and the right Riemann-Liouville fractional derivatives (RLFD) of order a of / 
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are defined, respectively, by 
and 

(1,14) 

Let AC([a,6]) represent the space of absolutely continuous functions on [a,b]. For n G N 
we denote by y4C"'[a, fe] the space of functions / which have continuous derivatives up to 
order ri — 1 on [a, b] such that /("-"i) g ACfa, b]. In particular, 74C"'^[a, 6] = AC[a, b]. 

The /e/^ and the right Caputo fractional derivatives (CFD) of order a G M.q of / G 
AC"'{[a,b]) are defined, respectively, by: 

Jn 

^D:/(x):=„jr"— /(x), x>a (1.15) 

and 

:= (-iy\J^--—f(x), X < b, (1.16) 

where n = [a] + 1. 

Remark 5. The Riemann-Liouville approach requires the initial conditions for differential 
equations in terms of non-integer derivatives, which are very difficult to be interpreted 
from the physical point of view [90], whereas the Caputo approach uses integer-order 
initial conditions that are more easy to get in real world problems (SH]. Caputo approach 
it's also preferable when we need fractional derivatives of constants to be zero. The Caputo 
fractional derivative is zero for a constant while the Riemann-Liouville derivative is not. 

Remark 6. If /(a) = /'(a) = ■ ■ ■ = /'■""^^(a) = 0, then both Riemann-Liouville and Caputo 
derivatives coincide. In particular, for a G (0, 1) and /(a) = one has '^D^f{t) = aD"/(t). 
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Chapter 2 
Time Scales 



In this chapter we recall some basic results on time scales (see Section |2TT]) that we use 
in the sequel. In Section 12.21 we review some results of the calculus of variations on time 
scales. 

2.1 Basic definitions 

Definition 7. A time scale is an arbitrary nonempty closed subset of M and is denoted 
by T. 

Example 8. Here we just give some examples of sets that are time scales and others that 
are not. 

• The set M is a time scale; 

• The set Z is a time scale; 

• The Cantor set is a time scale; 

• The set C is not a time scale; 

• The set Q is not a time scale. 

In applications, the quantum calculus that has, as base, the set of powers of a given 
number q is, in addition to classical continuous and purely discrete calculus, one of the most 
important time scales. Applications of this calculus appear, for example, in physics |70j 
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and economics |82]. The quantum derivative was introduced by Leonhard Euler and a 
fractional formulation of that derivative can be found in |96[[97]. 

In this thesis attention will be given to time scales Z, M and {h'L)a = {a, a + h,a + 
2/i, ...}, a G M, h>0. 

A time scale of the form of a union of disjoint closed real intervals constitutes a good 
background for the study of population (of plants, insects, etc.) models. Such models 
appear, for example, when a plant population exhibits exponential growth during the 
months of Spring and Summer, and at the beginning of Autumn all plants die while the 
seeds remain in the ground. Similar examples concerning insect populations, where all the 
adults die before the babies are born can be found in |1T |I111 . 



The following operators of time scales theory are used, in literature and throughout this 
thesis, several times: 

Definition 9. The mapping (j : T — )• T, defined by a{t) = inf {s G T : s > t} with inf = 
supT (i.e., a{M) = M if T has a maximum M) is called the forward jump operator. 
Accordingly, we define the backward jump operator p : T — )■ T by p{t) = sup {s G T : s < t} 
with sup0 = inf T (i.e., p(m) = m if T has a minimum m). The symbol denotes the 
empty set. 

The following classification of points is used within the theory: a point t G T is called 
right-dense, right- scattered, left-dense or left- scattered if a{t) = t, a{t) > t, p{t) = t, p{t) < 
t, respectively. A point t is called isolated if p(t) < t < a(t) and dense if p(t) = t = cr(t). 

Definition 10. A function / : T — )■ M is called regulated provided its right-sided limits 
exist (finite) at all right-dense points in T and its left-sided limits exist (finite) at all 
left-dense points in T. 

Example 11. Let T = {0} U |- : n G n| U {2} U |2 - - : G n| and define / : T ^ 
R by 

' if t 7^ 2; 



t ift = 2. 



m ■■= 

Function / is regulated. 

Now, let us define the sets T''", inductively: 

rjpK^ _ _ ^ . ^ non- maximal or left-dense} 
and T^" = ('T«"-^)«, n > 2. 
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Definition 12. The forward graininess function : T — )■ [0, oo) is defined by /^(t) = 
cr(t) - 1. 

Remark 13. Tlirougliout tlie tliesis, we refer to tlie forward graininess function simply as 
the graininess function. 

Definition 14. The backward graininess function z/ : T — )■ [0, oo) is defined by = 
t-p{t). 



Example 15. If T 



then a{t) = p{t) = t and = i/(t) = 0. If T = [-4, 1] UN, then 
a{t) = 



t ate [-4, 1) 

t + 1 otherwise. 



while 



Moreover, 



and 



Pit) 



flit) 



u{t) 



t iftG[-4,l]; 
t — 1 otherwise. 

if t e [-4,1); 

1 otherwise, 

if t G [-4,1]; 

1 otherwise. 



Example 16. If T = Z, then a{t) = t + 1, p{t) = t - 1 and p{t) = v{t) = 1. If T = (/iZ)„ 
then a{t) = t + h, p{t) = t — h and p{t) = u{t) = h. 

For two points a, 6 G T, the time scale interval is defined by 

[a,b]j = {t eT : a <t < b}. 

Throughout the thesis we will frequently write /'^(t) = f{a{t)) and /^(t) = f{p{t)). 
Next results are related with differentiation on time scales. 

Definition 17. We say that a function / : T — R is A-differentiable at t G T'' if there is 
a number f^it) such that for all e > there exists a neighborhood U oi t such that 



\nt)~fis)~f^it)ia{t)-s)\<e\ait) 



for all s G t/. 



We call f^{t) the ^-derivative of / at t. 
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The A-derivative of order n G N of a function / is defined by (t) = ( (t) j , 

t G T"", provided the right-hand side of the equahty exists, where = /. 
Some basic properties will now be given for the A-derivative. 

Theorem 18. Theorem 1.16] Assume / : T — )■ M a function and let t eT'^. Then 

we have the following: 

1. If f is A-differentiable at t, then f is continuous at t. 

2. If f is continuous at t and t is right-scattered, then f is A-differentiable at t with 

3. Ift is right-dense, then f is A- difjerentiable at t if and only if the limit 

s^t S — t 

exists as a finite number. In this case, 

^ ' s^t S-t 

4. If f is A-differentiable att, then 

r(t) = /(t)+/i(t)/^w. (2.2) 

It is an immediate consequence of Theorem [18] that if T = R, then the A-derivative 
becomes the classical one, i.e., = /' while if T = Z, the A-derivative reduces to the 
forward difference f^{t) = Af{t) = f{t + 1) - f{t). 

Theorem 19. [Jlj Theorem 1.20] Assume /, (7 : T — )■ M are A-differentiable at t E T'^ . 
Then: 

1. The sum / + ^ : T ^ R A-differentiable at t and (/ + g)^{t) = f^{t) + g^{t). 

2. For any number ^ G R, ^/ : T ^ R «s A-differentiable at t and (^/)^(t) = ^f^{t). 

3. The product fg:T^Wis A-differentiable at t and 

{fgtit) = f^mt) + nt)g\t) = f{t)g\t) + (2.3) 
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Definition 20. A function / : T — )■ R is called rd- continuous if it is continuous at right- 
dense points and if the left-sided limit exists at left-dense points. 

Remark 21. We denote the set of all rd-continuous functions by Crd or CrdCT), and the set 
of all A-differentiable functions with rd-continuous derivative by C]^ or Cjjj(T). 

Example 22. Let T = No U 1 1 - - : n e N I and define / : T ^ M by 



f{t) ■■- 



if t G N; 
t otherwise. 



It is easy to verify that / is continuous at the isolated points. So, the points that need our 
careful attention are the right-scattered point and the left-dense point 1. The right-sided 
limit of / at exist and is finite (is equal to 0). The left-sided limit of / at 1 exist and is 
finite (is equal to 1). We conclude that / is rd-continuous in T. 

We consider now some results about integration on time scales. 

Definition 23. A function F : T — M is called an antiderivative of / : T** — > M provided 
F^{t) = f{t) for all t e T«. 

Theorem 24. JJJl Theorem 1.74] Every rd-continuous function has an antiderivative. 

Let / : T'^ :— M be a rd-continuous function and let F : T — )■ M be an antiderivative of 
/. Then, the A-integral is defined by J^^ f{t)At = F{b) — F{a) for all a, 6 G T. 

One can easily prove gH Theorem 1.79] that, when T = M then £ fit) At = £ fit)dt, 
being the right-hand side of the equality the usual Riemann integral, and when [a, 6] fl T 
contains only isolated points, then 



'fit) At 



5^/i(t)/(t) ifa<6; 

t=a 

if a = 6; (2.4) 

P{b) 

-J2fiit)fit), ifa>b. 



t=a 

Remark 25. When T = Z or T = hZ, equation (12. 4p holds. 
The A — integral also satisfies 

r.a(t) 



fir)AT = fiit)fit), ^GT^ (2.5) 
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Theorem 26. ^ Theorem 1.77] Let a,b,c eT, ^ eR and f,g e Crd(T^) . Then, 



1- laifit) + aimt = \i fmt + g{t)At. 

3. jlf{t)At = -j^f{t)At. 

I j^fmt = i:fmt+i'fmt- 

5. j:fmt = o. 

6. If\f{t)\<g{t) on [a,h]^, then 



We now present the integration by parts formulas for the A-integral: 
Lemma 27. (cf. gl\ Theorem L77]) Ifa,beT and f\ge C].^, then 



Remark 28. For analogous results on V-integrals the reader can consult, e.g., |43| . 

Some more definitions and results must be presented since we will need them. We start 
defining functions that generalize polynomials functions in R to the times scales calculus. 
There are, at least, two ways of doing that: 

Definition 29. We define the polynomials on time scales., g^, /ifc : — > M for A; G Nq as 
follows: 





(2.6) 



(2.7) 



go{t,s) = ho{t,s 



) = 1 \/s,teT 
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Remark 30. If we let 



{t,s) denote, for each fixed s G T, the derivative of hk{t,s) with 



respect to t, then we have the following equalities: 

1. h^{t, s) = hk^i{t, s) for ken, te T^; 

2. gtit, s) = gk-iia{t), s) for keN, te 

3. gi{t, s) = hi(t, s) = t — s for all s,tET. 

Remark 31. It's not easy to find the explicit form of and for a generic time scale. 

To give the reader an idea what is the formula of such polynomials for particular time 
scales, we give the following two examples: 

Example 32. For T = M and /c G Nq we have: 



Example 33. Consider T = = {g^ : k e Z] U {0}. For k eNo and g > 1 we have: 



2.2 Calculus of Variations 

It is human nature to optimize. Typically, we try to maximize profit, minimize cost, 
travel to a destination following the smallest route or in the quickest time. The human 
propensity to optimize and the necessity to develop systematic tools for that has created a 
area of mathematics called optimization. Our main interest in this chapter is to introduce 
some concepts necessary to solve problems which involve finding extrema of an integral 
(functionals). The area of mathematics that deals with these problems is a branch of 
optimization called calculus of variations. Because this kind of calculus involve functionals 
it was called, earlier, functional calculus. Originally the name "calculus of variations" came 
from representing a perturbed curve using a Taylor polynomial plus some other term which 
was called the variation. Let us consider, for now, the following basic (but fundamental) 
problem: seek a function y E [a, b] such that 






fc-i 



(2.9) 




(2.10) 
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with a, b, Uai y;, G M and u, v) satisfying some smoothness properties. 

Remark 34. Although we write (I2.10p as a minimization problem, we can formulate it as 
a maximization problem using the fact to minimize is the same as to maximize 



We start by giving a simple application of the calculus of variations on the real setting 
and then we refer to some results of the calculus of variations on time scales, which includes 
as special cases the classical calculus of variations (T = M) and the discrete calculus of 
variations (T = Z). 

Example 35. Given two distinct points A = {a,yi) and B = 1/2) in the plane our 
task is to find the curve of shortest length connecting them. In childhood, all of us learn 
that the shortest path between two points is a straight line. That can be proved using, for 
example, the theory of calculus of variations. The formulation of this problem, in the sense 
of the calculus of variations, is to find the function y{-) that solves the following problem: 



In [37] Martin Bohner initiated the theory of calculus of variations on time scales. The 
rest of the section will be dedicated to present some results that will be necessary during 
the thesis. For more on the subject we refer to [51] and references therein. 

Definition 36. A function / defined on [a, b]j x M is called continuous in the second 
variable, uniformly in the first variable, if for each s > there exists 6 > such that 
\xi — X2I < 5 implies \ f{t,xi) — f{t,X2)\ < e for all t G [a, 6]t. 

Lemma 37 (cf. Lemma 2.2 in [37]). Suppose that F{x) = J^^ f{t,x)At is well defined. If 
fx is continuous in x, uniformly in t, then F'{x) = J^^ fx{t,x)At. 

Let us now extend (12.101) to a generic time-scale. 
Let a,b eT and u, v) : [a, b]j x — )• M. Find a function y E C]^ such that 




(2.11) 




(2.12) 



with ya, yt G M. 



Remark 38. If we fix T = M problem (|232|) reduces to flZTOj) . 
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Remark 39. If we fix T = Z we obtain the discrete version of fl2.12p . The goal is to find a 
function y E Cj^ such that 



6-1 



C[y{-)] = J2 L{t, y{t + 1), 2/^(t))At ^ min, y{a) = ya, y{h) = y^, (2.13) 



t=a 



where a, b E Z with a < b; ya, ?/b G M and L : Z x — M. 
Definition 40. For / G C]^ we define the norm 

11/11= sup \nt)\ + sup \f^it)\. 

A function y E Cl^ with y(a) = ya and y(b) = yb is called a (weak) local minimizer for 
problem (I2.12p provided there exists 6 > such that C{y) < C{y) for all y E C].^ with 
y{a) = ya and y{b) = yb and \\y - y\\ < S. 

Definition 41. A function r] E Cl^ is called an admissible variation provided r] ^ and 
77(a) = 77(6) = 0. 

Lemma 42 (cf. Lemma 3.4 in jS^). Let y^rj E C\.^ be arbitrary fixed functions. Put 
f{t,e) = L{t,y-{t) + £r/-(t),y^(t) + er]^{t)) and $(£) = C{y + er]), e E M. // /, and f,, 
are continuous in e, uniformly in t, then 

$'(5) = flL^it, y^{t), y'-itMit) + L,(t, y^{t), y\t))rj^{t)]At, 

$"(5) = f{L^M{tM{t)f + 2r/-(t)L„„[y](t)r/^(t) + L™M(i)(r/^(t))2}At, 

J a 

where [y]{t) = {t , y" {t) , y^ {t)) . 

The next lemma is the time scales version of the classical Dubois-Reymond lemma. 
Lemma 43 (cf. Lemma 4.1 in [33). Let g E Crd{[a,b]j) . Then, 

cb 

g{t)7]^{t)At = 0, for all t] E Cl^{[a, b]j) with r/(a) = r]{b) = 0, 
holds if and only if 

g(t) = c, on [a, b]^ for some c G M. 

The following theorem give first and second order necessary optimality conditions for 
problem fl232|) . 
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Theorem 44. Suppose that L satisfies the assumption of Lemma If V ^ C'^d « 

(weak) local minimizer for problem given by fl2.12p . then necessarily 

L L^[y](t) = Lu[y]{t), t G [a, (time scales Euler-Lagrange equation). 

2. L,,[y]{t) + ^^{t){2LuMit) + Kt)Luu[y]{t) + > 0, t e [a,fe]^' 

(time scales Legendre's condition), 

where [y]{t) = {t.y" {t),y^{t)) and a* = ^ if a e M\{0} and 0* = 0. 

Using the theory of calculus of variations on time scales and considering yi = and 
?/2 = 1 in Example [351 can generalize that example to the following one. 

Example 45. [STJ Example 4.3] Find the solution of the problem 

f y^l + (y^fAt min, y{a) = 0, y{b) = 1. (2.14) 
Writing fl2.14p in form of fl2.12p we have 



L{t, u,v) = Vl + Luit,u,v) = and Ly{t,u,v) 



1 + 

Suppose y is a local minimizer of fl2.14p . Then, by item 1 of Theorem HU equation 

L^{t,rit),y^{t)) = 0, tG[a,6]^ (2.15) 
must hold. The last equation implies that there exist a constant c G M such that 

L,{t,r{t),y^it)) = c, te[a,b], 



I.e., 



y^{t) = c^Jl + {y^)\ te[a,b] (2.16) 

holds. Solving equation (I2.16P we obtain 

y{t) = for all t G [a, 6], 

which is the equation of the straight line connecting the given points. 

We refer the reader to the PhD thesis [5T] for more recent results on calculus of vari- 
ations on time scales. Here our main interest on the subject is just as the starting point 
to go further and start developing the fractional case on the time scale Z (Chapter [3]), 
T = {h7j)^ (Chapter S]) and possible on an arbitrary time scale (cf. Chapter [S]). Before the 
original work of this thesis, already published in the international journals [33 | I34 ) . results 
on the fractional calculus of variations were restricted to the continuous case T = M (see, 

e.g., [HESllHo]). 
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Chapter 3 



Fractional Variational Problems in 



T = Z 



In this chapter we introduce a discrete-time fractional calculus of variations on the 
time scale Z. First and second order necessary optimality conditions are established. We 
finish the chapter with some examples illustrating the use of the new Euler-Lagrange and 
Legendre type conditions. 

3.1 Introduction 

In the last decades, considerable research has been done in fractional calculus. This is 
particularly true in the area of the calculus of variations, which is being subject to intense 
investigations during the last few years [25 | l26 [ l33 |ll02[ll03] . 

As mentioned in Section 11.11 Miller and Ross define the fractional sum of order u > 



Remark 46. This was done in analogy with the left Riemann-Liouville fractional integral 



by 




t-u 



(3.1) 



S = ' 



a 



of order u>0 (cf. ffTTU]) ) 
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which can be obtained via the solution of a hnear differential equation [86|[87]. Some basic 
properties of the sum in (13.11) were obtained in [86]. 

More recently, F. Atici and P. Eloe [201I2T] defined the fractional difference of order 
a > 0, i.e., A"/(t) = A™(A~(™"°')/(t)) with m the least integer satisfying m > a, and de- 
veloped some of its properties that allow to obtain solutions of certain fractional difference 
equations. 

To the best of the author's knowledge, the first paper with a theory for the discrete 
calculus of variations (using forward difference operator) was written by Tomlinson Fort |54) 
in 1937. Some important results of discrete calculus of variations are summarized in [7T| 
Chap. 8]. It is well known that discrete analogues of differential equations can be very 
useful in applications pB | [7T ] . Therefore, we consider pertinent to start here a fractional 
discrete-time theory of the calculus of variations. 

Throughout the chapter we proceed to develop the theory of fractional difference cal- 
culus, namely, we introduce the concept of left and right fractional sum/difference (cf. 
Definitions HT] and |53] below) and prove some new results related to them. 

We begin, in Section 13. 2^ to give the definitions and results needed throughout. In 
Section 13.31 we present and prove the new results; in Section 13.41 we give some examples. 
Finally, in Section 13.51 we mention the main conclusions of the chapter, and some possible 
extensions and open questions. Computer code done in the Computer Algebra System 
Maxima is given in Appendix A. 

3.2 Preliminaries 

We begin by introducing some notation used throughout. Let a be an arbitrary real 
number and b = k + a foT a certain /c G N with k > 2. In this chapter we restrict ourselves 
to the time scale T = {a, a + 1, . . . ,b}. 

Denote by J-" the set of all real valued functions defined on T. According with Exam- 
ple Uni we have a(t) = t + 1, p{t) = t — 1. 

The usual conventions Yl't=c /(^) = 0, c G T, and Y[7=o /(^) ~ ^ remain valid here. As 
usual, the forward difference is defined by Af{t) = f^{t) — f{t). If we have a function 
/ of two variables, f{t,s), its partial (difference) derivatives are denoted by At and A^, 
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respectively. Recalling (11. 8p we can, for arbitrary x,y eM, define (when it makes sense) 

x(^) = + , (3.2) 

r(x + l-y) ^ ^ 

where F is the gamma function. 

While reaching the proof of Theorem I^S] we actually "find" the definition of left and 
right fractional sum: 

Definition 47. Let / G J-". The left fractional sum and the right fractional sum of order 
z/ > are defined, respectively, as 

^ t—i/ 

a^rfit) = - ^(^))^^"'V(^), (3.3) 

^ ^ s=a 

and 

tAr/(t) = E - ^{t)Y'-'^f{s). (3.4) 

Remark 48. The above sums (13. 3p and (13. 4p are defined for t E {a + u, a + u + 1 , . . . , b + u} 
and t G {a — u, a — u + 1, . . . , b—u}, respectively, while f{t) is defined for t G {a, a+1, . . . , 6}. 
Throughout we will write (13. 3p and (13. 4p . respectively, in the following way: 

Ar/(t) = + ^ - ^(^))^'^-^V(^), t G T, 



Remark 49. The left fractional sum defined in (13.31) coincides with the fractional sum 
defined in [86] (see also (II. 9p ). The analogy of (13. 3 p and (13. 4 p with the Riemann-Liouville 
left and right fractional integrals of order z/ > is clear: 

al-'fix) = ^ J\x - Sy-'f{s)ds, 

./r/(x) = ^ l\s-xr'fis)ds. 



25 



CHAPTER 3. FRACTIONAL VARIATIONAL PROBLEMS IN T = Z 



It was proved in [86] that limj,_j.o '^/(t) = f{t). We do the same for the right 
fractional sum using a different method. Let z/ > be arbitrary. Then, 

b 

^ ' s=t 

^ ^ s=a(t) 

S = (j{t) 

s=a{t) ^ ' 

Therefore, \\m.y^{) iJ\^^ ] if) = f{t)- It is now natural to define 

^A'^J{t)=A'j{t) = f{t), (3.5) 

which we do here, and to write 

,Ai'^f(t) = f{t) + —^J2(t + u-a{s)Y''-'^f{s), teT, u>0, (3.6) 

^ ' s=a 

b 

*Ar/(i) = fit) + YrTT^ E + ^ - ^it)Y'-'^fi^), ^ e T, u>o. 

' s=a{t) 



The next theorem was proved in [21} . 
Theorem 50. /H]/ Let / G J-" and u > 0. Then, the equality 



A/(t) = A(„Ar/(t)) - ^^^^^Y(^^f(^)^ ^ ^ 



holds. 



Remark 51. It is easy to include the case = in Theorem [501 Indeed, in view of (11. 6p 
and (13. 5p . we get 



A^Afit) = A(.A-V(t)) - fT^At + ^ - a)^^~'^f{a), t G T^ (3.7) 



for all u > 0. 



26 



3.2. PRELIMINARIES 



Now, we prove the counterpart of Theorem [50] for the right fractional sum. 
Theorem 52. Let / G J-" and u > 0. Then, the equality 

V 

r(^ + i) 

holds. 



A;(^)A/(t) = ^,r^,{h + V - a{t)t-'^]{h) + t G T^ (3.8) 



Proof. We only prove the case z/ > as the case z/ = is trivial (see Remark I^T]) . We start 
by fixing an arbitrary t G T*^ and prove that for all s G T'', 

A, {{s + u- a(t))('^-i)/(.)) = (zy - l)(s + z/ - a{t)Y''-'^ris) + {s + u- a{t)Y''-'^Af{s). 

(3.9) 

By definition of forward difference with respect to variable s we can write that 

A, {{s + u- a{t)Y''~''>f{s)) = {a{s) + v - cr{t)f-'^ T (s) - {s + u - a{t)f-'^ f{s) 

= {s + u- tt'^^ ris) -{s + u-t- f{s) 
r(s + z/-t + l)r(s) T{s + u-t)f{s) 



r(s-t + 2) r(s-t + l) 



(3.10) 



Applying (11.61) to the numerator and denominator of first fraction we have that fl3.10p is 
equal to 

{s + u~t)r{s + u-t)r{s) _ r{s + u-t)f{s) 
(s-t + i)r(s-t + 1) r(s-t + i) 

[(g + t)r{s) -{s-t + l)f{s)] T{s + u- t) 
{s-t + l)T{s-t + l) 

[{s-t + l)r (s) + (z/ - l)r (s) -{s-t + l)f{s)] T{s + u- t) 
{s-t+l)T{s-t + l) 

T[s + u — a[t) — (z/ — 2) + 1) T[s + u — a[t) — (z/ — 1) + 1) 

= (z. - l)(s + z. - a{t)t-^^r{.s) + {s + u- a(t))(^-i)A/(.) . 
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Now, with fl3.9p proven, we can state that 



s=t 



fis) 



s=b ^ b-l 



5^(z.-l)(. + z.-a(t))(-2)r(.) 



b-l 



We now compute A((A^''/(t)): 



r(^) 



s=a{t) 



1) 



Til.) 



s=t 

b 

s=a{t) 

b 

b 

s=a{t) 



i)(- 



-fit) 



fit) 



Tiu) 



Since t is arbitrary, the theorem is proved. 



□ 



Definition 53. Let < a < 1 and set yU = 1 — a. Then, the left fractional difference and 
the right fractional difference of order a of a function f ^ T are defined, respectively, by 



and 



,A^/(t) = -AGA//(t))), t G T^ 
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3.3 Main results 

Our aim is to introduce the discrete-time (in time scale T = {a, a + 1, . . . , &} ) fractional 
problem of the calculus of variations and to prove corresponding necessary optimality con- 
ditions. In order to obtain an analogue of the Euler-Lagrange equation (cf. Theorem [59]) 
we first prove a fractional formula of summation by parts. Our results give discrete ana- 
logues to the fractional Riemann-Liouville results available in the literature: Theorem [5^ is 
the discrete analog of fractional integration by parts |1051I108] : Theorem [591 is the discrete 
analog of the fractional Euler-Lagrange equation of Agrawal [3l Theorem 1]; the natural 
boundary conditions fl3.24p and ( I3.25P are the discrete fractional analogues of the transver- 
sality conditions in [Hllin]- However, to the best of the author's knowledge, no counterpart 
to our Theorem |62] exists in the literature of continuous fractional variational problems. 

3.3.1 Fractional summation by parts 

The next lemma is used in the proof of Theorem [551 

Lemma 54. Let f and h be two functions defined on and g a function defined on 
X T". Then, the equality 

6-1 T-1 6-2 6-1 

Y.f{r)Y,9{r,s)h{s) = Y,h{r) 9{s,r)f{s) 

T=a s=a T=a s=u{t) 

holds. 

Proof. Choose T = Z and -F('r, s) = f{T)g{T, s)h{s) in Theorem 10 of |9]. □ 

The next result gives a fractional summation by parts formula. 

Theorem 55 (Fractional summation by parts). Let f and g be real valued functions defined 
on and T, respectively. Fix < a < 1 and put n = 1 — a. Then, 



6-1 6-2 

f{t)aA1git) = fib - l)gib) - fia)g{a) + ,A^(,)/(t)^7'^(t) 

t=a t=a 

( 6-1 6-1 
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Proof. From f)3.7|) we can write 



J2fit)aAtg{t) = J2fmUA;^9{t)) 



t=a 



t=a 
6-1 



t=a 
b-l 



r(/i + i) 



(t + /i-a)(^-^)^(a) 



J2 maA;^Agit) + J2 f77^(^ + - (3.11) 



t=a 



t=a 



r(/x + i) 



Using (13. 6p we get 



6-1 



t=a 



6-1 



6-1 i-1 



6-1 



r(/i + i) 
/i 



t=a s=a 
6-2 6-1 



t=a ' t=a s=a-{t) 



t=a 
6-2 



f{h - l)[g{h) - g{b - 1)] + E Ag{t),A;i:)f{t), 



t=a 



where the third equahty follows by Lemma [Ml We proceed to develop the right hand side 
of the last equality as follows: 



6-2 



fib - 1Mb) - gib - 1)] + E A^?(t),A;(y (t) 



t=a 



fib-i)[gib)-gib-i)]+ ^(t),A-y(t)|^^^ -Y,fm{tA;i;^fit)) 



t=b-l 



6-2 



t=a 



fib-l)gib)-fia)gia) 



fx 



6-1 



r(/i + i) 



gia) E is + f^-aia)Y^-'^fis) 



s=a(a) 



6-2 



t=a 



where the first equality follows from the usual summation by parts formula. Putting this 
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into (13. lip , we get 



6-2 



E m^^tgit) = fib - l)gib) - f{a)g{a) + (*A°(6)/(^)) g^it) 



t=a 



t=a 



g{a)fi ^ (t + /x-a)(^-i) 



+ 



r(/i + 



X — ^ 



t=a 



r(/i) 



fit) 



gia)f^ 
r(/i + i) 



6-1 



J2 is + f^-aia)Y^-'^fis). 



s=cr{a) 



The theorem is proved. 



□ 



3.3.2 Necessary optimality conditions 

We begin to fix two arbitrary real numbers a and /3 such that a, /3 G (0, 1]. Further, 
we put fi = 1 — a and z/ = 1 — /3. 

Let a function L{t, u,v,w) : T'' x M x M x M — )■ M be given. We assume that the second- 
order partial derivatives Luu, Luv, Luw, L^^, L„„, and Ly^^ exist and are continuous. 

Consider the functional £ : J-" — t- M defined by 

6-1 

Civi-)) = J2Lit,y-{t),,Atyit),Xyit)) (3-12) 

t=a 

and the problem, that we denote by (P), of minimizing (I3.12p subject to the boundary 
conditions y{a) = A and y{b) = B (A, 5 G M). Our aim is to derive necessary conditions 
of first and second order for problem (P). 

Definition 56. For / G J-" we define the norm 

11/11 = max \nt)\ + max UAr/(t)| + max |,Af/(t)|. 

A function y ^ T with y[a) = A and y{b) = B is called a local minimizer for problem 
(P) provided there exists d > such that Ciy) < Ciy) for all ?/ G J-" with y{a) = A and 
yib) = B and \\y - y\\ < 5. 

Remark 57. It is easy to see that Definition [35] gives a norm in J-". Indeed, it is clear that 
11/11 is nonnegative, and for an arbitrary / G J-" and A; G M we have ||A;/|| = |A;|||/||. The 
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triangle inequality is also easy to prove: 

11/ + g\\ = max + git)\ + max UA°(/ + g){t)\ + max |,A°(/ + gm 

iGT*^ teT« teT" 

< max[|/(t)| + \g{t)\] + max[UAtf{t)\ + 
+ max[|,A^/(t)| + |,A^^7WI] 



< 11/11 + 11^^11- 

The only possible doubt is to prove that ||/|| =0 implies that f{t) = for any t E T = 
{a,a + 1, . . . ,b}. Suppose 1 1 /I I = 0. It follows that 

max|r(t)|=0, (3.13) 

maxUA,"/(t)|=0, (3.14) 

max|,Af/(t)| = 0. (3.15) 

From (I3.13P we conclude that /(t) = for all t G {a + 1, . . . , 6}. It remains to prove that 
/(a) = 0. To prove this we use fl3.14p (or f l3.15p ). Indeed, from f l3.13p we can write 

a^?f{t) = A + 1 - " - ^(^))^""V(^) 

, /m t 

5^(t + 2 - a - a(s))(-")/(.) - ^(t + 1 - a - a(s))(-")/(s) 

\s=a s=a 

{{t + 2-a- a(a))("°V(a) - {t + 1 - a - a(a))("")/(a)) 



r(l - a) 
1 



r(l - a) 
•^(«) -A(t + l-a-a(a))("") 



r(l - a) 

and since by fl3.14p aA"/(t) = 0, one concludes that /(a) = (because (t + 1 — a — a (a) )'•""•' 
is not a constant). 

Definition 58. A function 77 G J-" is called an admissible variation for problem (P) provided 
rj ^ and 77(a) = 77(6) = 0. 

The next theorem presents a first order necessary optimality condition for problem (P). 

Theorem 59 (The fractional discrete-time Euler-Lagrange equation). If y E T is a local 
minimizer for problem (P), then 

L^it) + iA^(,)L4y](t) + aAf L^[y](t) = (3.16) 
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holds for all t G , where the operator [■] is defined by 



[y]{s) = {s,y'^{s),^A':y{s),Xyis)). 



Proof. Suppose that y{-) is a local minimizer of £(■). Let ?7(-) be an arbitrary fixed admis- 
sible variation and define the function $ : 



M-)V M-)\ 



^ M by 



^e) = C{y{-)+evi-)). 



(3.17) 



This function has a minimum at e = 0, so we must have $'(0) = 0, i.e. 

b-l 

J2 [LuimKit) + L4y]{t)aKv{t) + L^mt^hit) 

t=a 

which we may write, equivalently, as 



b-2 



t=a 



6-1 



t=a 



t=a 



Using Theorem [55| and the fact that 77(a) = 77(6) = 0, we get for the third term in (I3.18P 
that 



fe-i 



b-2 



t=a 



t=a 

Using (13. 8p it follows that 

6-1 

6-1 

-5^L4y](t)A(,Arr/(t)) 



(3.19) 



t=a 



t=a 
6-1 



t=a 
'6-1 



A,T,)Ar/(t) 



uri{b) 



(3.20) 



6-1 



, t=a 



t=a 
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We now use Lemma 15^ to get 



6-1 



^L4y](t),A;(^)Ar/(t) 



t=a 



6-1 6-2 6-1 

V 



6-1 6-1 t-1 

V 



t=a ^ ' t=a s=a 

6-1 

t=a 

We apply again the usual summation by parts formula, this time to fl3.2ip . to obtain: 

6-1 

5^Ar/(t),A-^L^[y](t) 

t=a 

6-2 

= ^ Ar/(t),A7'^L^[y](t) + (r/(6) - r/(p(6))),A7'^L„,[y](t)|,=,(,) 

t=a 

t=a 

+ rj{b)aA-''L^[y]{t)\t=pib) - v{b - l)aAi'' L^[y]{t)\t=p(b) 

6-2 

= r/(6),A-^L,4y](t)|i=,(,) - r]{a)aA;''L4y]{t)\t=a - 5^r/'^(t)aAfL^[y](t). 

t=a 

Since r/(a) = = it follows, from (IS:^ and (IX^ . that 

6-1 6-2 

and, after inserting in fl3.20p . that 

6-1 6-2 

YL^ymtA^bVit) = 5^r/-(t).Af^4^]W- (3.23) 

t=a t=a 

By flCTjl and flX^ we may write flXTH]) as 

6-2 

Y [Lu[m) + *A^(,)L4y](t) + .AfL,^[y](t)] = 0. 

t=a 

Since the values of r]'^{t) are arbitrary for t G T''^, the Euler-Lagrange equation f l3.16p 
holds along y. □ 
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Remark 60. If the initial condition y{a) = A is not present (i.e., y{a) is free), we can use 
standard techniques to show that the following supplementary condition must be fulfilled: 



'6-1 



- Ua) + ^^^^^ lEit + f,- aY'^-'^LMit) 



, t=a 

b-1 



^ (t + ^ _ a{a)Y^-'^LM{t) I + LM = 0. (3.24) 

t=a{a) 



Similarly, if y{b) = B is not present (i.e., y{b) is free), the equality 



. t=a 

b-2 



- E(/^(^) + " - ^{t)Y''''^LMit) = (3.25) 



t=a 



holds. We just note that the first term in (13.251) arises from the first term on the left hand 
side of fl3.18p . Equalities fl3.24p and f l3.25p are the fractional discrete-time natural boundary 
conditions. 

The next result is a particular case of our Theorem [591 

Corollary 61 (The discrete-time Euler-Lagrange equation - cf., e.g., [37|I52|). If y is a 

solution to the problem 

6-1 

C{y{-)) = y^L{t,y{t + l),^y{t))-^min 

(3.26) 

y{a) = A, y{b) = B, 
then Luit, y{t + 1), Ay(t)) - AL^(t, y{t + 1), Ay(t)) = Q for all t e {a, . . . ,b - 2} . 
Proof. Follows from Theorem El] with a = 1 and a L not depending on w. □ 

We derive now the second order necessary condition for problem (P), i.e., we obtain 
Legendre's necessary condition for the fractional difference setting. 
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Theorem 62 (The fractional discrete-time Legendre condition). If y E is a local mini- 
mizer for problem (P), then the inequality 

Luu[y\{t) + 2L^,[y\it) + L,M(t) + - 1)' 



+ 2(1. - l)L^\y\(t) + 2{p - l)L„Jg]{a(t)) + L^4y]{t)(l - vf 



holds for allte where [y]{t) = (t, rW, aA°y(t), tAfy(t)). 

Proof. By the hypothesis of the theorem, and letting $ be as in f l3.17p . we get 

$"(0) > (3.27) 
for an arbitrary admissible variation ?7(-). Inequality (I3.27P is equivalent to 

[L^^miv^it))' + 2L^S]{t)v''{t),A1r^{t) + L,4y](t)(„A°r^(t))^ 



t=a 



+2L^4ymv''{t)tA'bV{t) + 2L,4y]{t)aAtv{t)tA'^vit) + L^4y]{t){tA^,r]{t)f 



> 0. 



2 

Let r G T'^ be arbitrary and define 77 : T — )■ M by 



1 if t = (T(r); 
otherwise. 



It follows that 77(a) = 77(6) = 0, i.e., rj is an admissible variation. Using (13. 7p (note that 
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ri{a) = 0), we get 



6-1 



t=a 



h-i 



t=a 



r(/i + i) 



Luu[y]{r) + 2Luv[y]{r) + L^„[y](r) 

/ t-l ^2 

Ar^(t) + V(t + /i - a(.))(^-i)Ar/(s) 



t=a(T) 



Observe that 



i=<7{a(r)) 



6-1 / t-l 



t=a(T) 



We show next that 



6-1 



t-l 



t=<7(<7(r)) 



r(/i + i) 



5^(t + /i-a(.))(^"^)Ary(s) 



6-1 



t=<7(f7(r)) 



_ 7-)r(t - r) 
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Let t e [a{a{T)), 6 - 1] n Z. Then, 



t-i 



r(/i + i) 
r(/i + i) 
r(/^ + i) 



s=a s=(7(r) 

[(t + /i - a(r))(^-i) - (t + /i - a(a(r)))(^-i)] 

r(t + ^-(7(r) + 1) r(t + /i-a(a(r)) + 1) 



r(t + /X - a(r) + 1 - (/X - 1)) r(t + /X - (T(a(r)) + 1 - (/x - 1)) 



r(/^ + i) 
r(/i + i) 



r(t + /i-r) _ r(t-r + ;u - 1) 

r(t-r + i)~ r(t-r) 

(t + ;U - r - l)r(t + /i - r - 1) {t- r)r(t - r + /i - 1) 



(i-r)r(t-r) 



(3.28) 



{t-T)T{t-T) 
fi (/i — l)r(t — r + /i — 1) 
r(/i + l) (t - r)r(t - r) 

M/^-i)n-x'(/i+^ + i) 

(t - r)r(t - r) 



which proves our claim. Observe that we can write t^bVit) = ^t^p(b)^'7(^) since 77(6) = 0. 
It is not difficult to see that the following equality holds: 



6-1 



J22L^4y]itKit)tA^bV{t) = -5^2L,^[y](t)r/'^(t)tA;(^,)Ar/(t) = 2L„„,[y](r)(i. - 1). 



t=a 



t=a 



Moreover, 



b~l 



Y,2L,n.[y]it)aAtVit)t^Ut) 



t=a 



t=a 



r(^ + i) 



IX — (T{S 



))('^"i)Ar/(^ 



b-i 



Ar]{t) + 



=a{t) 



2{y - l)L,^[y](r) + 2(/i - l)L„^[y](a(r)). 
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Finally, we have that 



6-1 



t=a 



t=a 

T-1 



Ar]{t) 



V 



s=a(t) 



t=a 



V 



h-\ 



s=(y(t) 

+ L^^[y](r)(l - vf + L^^[y](o-(r)) 

t=a L I "T / 

[y]{T){l-uf + L^4y]{a{T)). 
Similarly as we have done in (13. 28 p . we obtain that 

V 



r 



(z/+ 1) L (^(^) _ t)r(a(r) - t) 



We are done with the proof. 



□ 



A trivial corollary of our result gives the discrete-time version of Legendre's necessary 
condition. 

Corollary 63 (The discrete-time Legendre condition - cf., e.g., [STflM]). Ij y is a solution 
to the problem (I3.26p . then 

Luu[y]{t) + 2L„,[y](t) + K,[y]{t) + L,,[y]{a{t)) > 

holds for all t e T"', where [y]{t) = {t,y''{t),Ay{t)). 

Proof. We consider problem (P) with a = 1 and L not depending on w. The choice a = 1 
implies /i = 0, and the result follows immediately from Theorem [62j □ 



3.4 Examples 

In this section we present three illustrative exarnples. The results were obtained using 
the open source Computer Algebra System Maximao All computations were done running 
^http : / /maxima. sourcef orge .netl 
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Maxima on an Intel® Core™2 Duo, CPU of 2.27GHz with 3Gb of RAM. Our Maxima 
definitions are given in Appendix lAl 

Example 64. Let us consider the following problem: 

6-1 

Ja(y) = 5Z(oA,"y(t))'-^min, y{0)=A, y{h) = B . (3.29) 
t=o 

In this case Theorem ES] is trivially satisfied. We obtain the solution y to our Euler- 
Lagrange equation ( I3.16P for the case b = 2 using the computer algebra system Maxima. 
Using our Maxima package (see the definition of the command extremal in Appendix 
we do 



Ll:v~2$ 

extremal (LI , , 2 , A , B , alpha , alpha) ; 



to obtain (2 seconds) 



2aB + {a^-a^ + 2a) A 
yii) = • (3.30) 



For the particular case a = 1 the equality fl3.30p gives y(l) = which coincides with 

the solution to the (non-fractional) discrete problem 

1 1 

J](Ay(t))2 = 5^(y(t + l)-y(t))^-^min, 1/(0) = A, yi2)=B. 
t=o t=o 

Similarly, we can obtain exact formulas of the extremal on bigger intervals (for bigger 
values of b). For example, the solution of problem (I3.29P with 6 = 3 is (35 seconds) 

, , (6a'^ + 6a) B + (2a^ + 2a^ + 10a^-2a'^ + 12a) A 
1/(1) = 



^(2) 



3a^ + 6a^ + 15a^ + 12 
12«3 + 12^2 + 24a) B + {a^ + + 7 - + Aa"^ + 12 a) A 
6a^ + 12a^ + 30 + 24 
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and the solution of problem fl3.29p with 6 = 4 is (72 seconds) 

3 + 15 + 57 + 69 + 156 -12a^ + 144 a 



where 



24 a^ + 72a^ + A8a^ 
+ z B. 



A 



a 



A 



5 + 22 + 32 + 67 + 35 + 54 + 72 a 
24 + 72 ^ 120 + 72 a ^ 

+ 1 ^- 

+ 6 a® + 30 + 60 + 117 + 150 - 4 + 216 + 288 



C 



^ = 4 + 24 + 88 + 120 + 196 a' 



144, 



a 



c 

72 + 288 + 792 + 576 + 864 a ^ 
+ : B. 



A 



UA a' 



528 + 720 + 1176 + 864 . 



C = 24a'^ 

Consider now problem fl3.29p with b = A, A = 0, and 5 = 1. In Table 13.11 we show the 
extremal values y{l), y{2), y(3), and corresponding Jq,, for some values of a. Our numerical 
results show that the fractional extremal converges to the classical (integer order) extremal 
when a tends to one. This is illustrated in Figure [3711 The numerical results from Table \37l\ 



and Figure [3^ show that for this problem the smallest value of Jq,, a g]0,1], occur for 
a = 1 (i.e., the smallest value of J„ occurs for the classical non-fractional case). 



a 


m 


m) 


m 




0.25 


0.10647146897355 


0.16857982587479 


0.2792657904952 


0.90855653524095 


0.50 


0.20997375328084 


0.35695538057743 


0.54068241469816 


0.67191601049869 


0.75 


0.25543605027861 


0.4702345471038 


0.69508876506414 


0.4246209666969 


1 


0.25 


0.5 


0.75 


0.25 



Table 3.1: The extremal values y{l), y{2) and y{3) of problem f l3.29p with b = 4, A = 0, 
and B = 1 for different a's. 



Example 65. In this example we generalize problem (I3.29P to 



2 2 

^,/3 = 5Z7i(oA,"y(t)) +72(tAfy(t)) ^min, y{0) = A, y{b) = B . (3.31) 
t=o 
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1.0 - 
0.9 
0.8 - 
0.7 
0.6 - 
0.5 + 
0.4 - 
0.3 
0.2 
0.1 
' 





+ 

X 



Figure 3.1: Extremal y{t) of Example [641 with b = 4, A = 0, B = 1, and different a's 
a = 0.25; x: a = 0.5; +: a = 0.75; *: a = 1). 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



Figure 3.2: Function of Example [64] with b = 4, A = 0, and B = 1. 
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As before, we solve the associated Euler-Lagrange equation fl3.16|) for the case 6 = 2 with 
the help of our Maxima package (35 seconds): 



L2 : (gamma [1] ) *v"2+ (gamma [2] ) *w'-2$ 
extremal (L2 , , 2 , A , B , alpha , beta) ; 

(2 72 /3 + 7i _ «2 ^ 2 7i a) A + (72 



72/32 + 272/3 + 271 «) 5 



272/32 + 27102 + 272 + 271 
Consider now problem (I3.3ip with 71 = 72 = 1, 6 = 2, A = 0, i? = 1, and /3 = a. In 
Table [3l2] we show the values of y{l) and Jq, := Ja^aivi^)) for some values of a. We con- 
cluded, numerically, that the fractional extremal y{l) tends to the classical (non- fractional) 
extremal when a tends to one. Differently from Example [Ml the smallest value of J^, 
a g]0, 1], does not occur here for a = 1 (see Figure 13. 3p . The smallest value of J^, 
a g]0, 1], occurs for a = 0.61747447161482. 



a 


m 


Ja 


0.25 


0.22426470588235 


0.96441291360294 


0.50 


0.375 


0.9140625 


0.75 


0.4575 


0.91720703125 


1 


0.5 


1 



Table 3.2: The extremal y{l) of problem fl3.3ip for different values of a (71 
6 = 2, A = 0, = 1, and /3 = a). 



72 
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Example 66. Our last example is a discrete version of the fractional continuous problem [6l 
Example 2]: 

J, = J]-(oA,"y(t))'-y-(t) -^min, 1/(0) = 0, 1/(2) = . (3.32) 



t=o 



The Euler-Lagrange extremal of f l3.32p is easily obtained with our Maxima package (4 
seconds) : 



L3: (l/2)*v"2-u;$ 

extremal (L3 ,0,2,0,0, alpha , beta) ; 



(3.33) 



For the particular case a = 1 the equality (13.331) gives yiV) = i, which coincides with the 
solution to the non-fractional discrete problem 

J2 ^ i^yit)? - y^it) = 5^ - {y{t + 1) - y{t)f - y{t + 1) min , 



t=o 



t=o 



1/(0) = 0, 1/(2) = 0. 



In Table [23] we show the values of y{l) and for some a's. As seen in Figure 
a = 1 one gets the maximum value of Jq, a g]0, 1]. 



for 



a 


m 


J a 


0.25 


0.94117647058824 


-0.47058823529412 


0.50 


0.8 


-0.4 


0.75 


0.64 


-0.32 


1 


0.5 


-0.25 



Table 3.3: Extremal values y{l) of (I3.32p for different a's 



3.5 Conclusion 

In this chapter we introduce the study of fractional discrete-time problems of the calcu- 
lus of variations of order a, < a < 1, with left and right discrete operators of Riemann- 
Liouville type. For a = 1 we obtain the classical discrete-time results of the calculus of 
variations |71] . 
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Main results of the chapter include a fractional summation by parts formula (Theo- 
rem [SS]), a fractional discrete-time Euler-Lagrange equation (Theorem [5^ . transversality 
conditions ( I3.24p and ( I3.25p . and a fractional discrete-time Legendre condition (Theo- 
rem [62]). From the analysis of the results obtained from computer experiments, we con- 
clude that when the value of a approaches one, the optimal value of the fractional discrete 
functional converges to the optimal value of the classical (non-fractional) discrete problem. 
On the other hand, the value of a for which the functional attains its minimum varies with 
the concrete problem under consideration. 



3.6 State of the Art 

The discrete-time calculus is also a very important tool in practical applications and 
in the modeling of real phenomena. Therefore, it is not a surprise that fractional discrete 
calculus is recently under strong development (see, e.g., 15 ^ 14 ^ 1 ^ 1^ and references 
therein). There are some recent theses [6711109] with results related with discrete calculus on 
time scale T = {a, a + l,a + 2, . . .}. Thesis |109] links the use of fractional difference equa- 



tions with tumor growth modeling, and addresses a very important question: by changing 
the order of the difference equations from integers to fractional, in which conditions, we 
are able to provide more accurate models for real world problems? For more on the subject 



see mm. 



In this chapter we use fractional difference operators of Riemann-Liouville type. However, 
there are some results related with fractional difference operators in Caputo sense [17] like 
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in [l5], where the existence of solutions for I VP involving nonlinear fractional difference 
equations is discussed. 

The results of this chapter are published in [23] and were presented by the author at 
the Workshop on Control, Nonsmooth Analysis and Optimization, Celebrating Francis 
Clarke's and Richard Vinter's 60th Birthday, Porto, Portugal, May 4-8, 2009 in a con- 
tributed talk entitled Necessary Optimality Conditions for Fractional Difference Problems 
of the Calculus of Variations. 
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Chapter 4 



Fractional Variational Problems in 

T = {hZ)a 

We introduce a discrete-time fractional calculus of variations on the time scale {h7j)a, 
h > 0. First and second order necessary optimality conditions are established. Examples 
illustrating the use of the new Euler-Lagrange and Legendre type conditions are given. 

4.1 Introduction 

Although the fractional Euler-Lagrange equations are obtained in a similar manner as 
in the standard variational calculus |102| . some classical results are extremely difficult to be 
proved in a fractional context. This explains, for example, why a fractional Legendre type 
condition is absent from the literature of fractional variational calculus. In this chapter we 
give a first result in this direction (cf. Theorem [M]) . 

Despite its importance in applications, less is known for discrete-time fractional sys- 
tems |102] . Our objective is proceed to develop the theory of fractional difference calculus, 
namely, we introduce the concept of left and right fractional sum/difference (cf. Defini- 
tion [TH])- In Section 112] we introduce notations, we give necessary definitions, and prove 
some preliminary results needed in the sequel. Main results of the paper appear in Sec- 
tion we prove a fractional formula of /i-summation by parts (Theorem [85]) . and neces- 
sary optimality conditions of first and second order (Theorems [88] and |9T| respectively) for 
the proposed /i-fractional problem of the calculus of variations fl4.15p . Before the end of 
the chapter we present in Section U3] some illustrative examples and give some conclusions 
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in Section 14.51 and the state of the art in Section 14.61 



4.2 Preliminaries 



One way to approach the Riemann-Liouville fractional calculus is through the theory 



of linear differential equations |100| . Miller and Ross j86] use an analogous methodology 
to introduce fractional discrete operators for the case T = Z = {a, a + l,a + 2, ...},aG]R. 
That was our starting point for the results presented in Chapter [31 Because the results 
obtained in the previous chapter were good, we start thinking how to generalize them to 
an arbitrary time scale. However, since we found several difficulties when we try to get 
results for a time-scale where the graininess function is not constant, here we keep the 
graininess function constant but not, necessarily, equal to one and we give a step further: 
we use the theory of time scales in order to introduce fractional discrete operators to the 
more general case T = {KL)a = {a, a + h,a + 2h, . . .}, a G M, /i > 0. 

Before going further we need to collected more results on the theory of time-scales. 

Until we state the opposite, the following results are true for a arbitrary time scale. 

For n G No and rd-continuous functions : T — ^ R, 1 < i < n, let us consider the nth 
order linear dynamic equation 



Ly = , where Ly = y 



A" 



n 

1=1 



A"- 



(4.1) 



A function y : T — t- M is said to be a solution of equation (14. ip on T provided y is n times 
delta differentiable on T'^" and satisfies Ly{t) = for all t G T'^". 



Lemma 67. gl[ p. 239] If z = [zi, 



Zn) -.T ^ W satisfies for all t G T'^ 
/ 1 ... \ 



= A{t)z{t), where A 



1 





V -Pn 





1 

-P2 -Pl ) 



(4.2) 



then y = z\ is a solution of equation (14. ip . Conversely, if y solves ( 14. ip on T, then 
z= (y,y^,..., y^"'') : T ^ R satisfies for all t G T''" 
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Definition 68. [HI p. 239] We say that equation (14. ip is regressive provided / + n{t)A{t) 
is invertible for all t G T'^, where A is the matrix in (14.21) . 

Definition 69. HH p. 250] We define the Cauchy function y : T x T^" ^ R for the 
linear dynamic equation (14. ip to be, for each fixed s G T''" , the solution of the initial value 
problem 



Tl-l 



Ly = 0, {a{s),s)=0, 0<t<n-2, " s) = 1 . (4.3) 

Theorem 70. JJl, p. 251] Suppose {yi, . . . , ?/,„} is a fundamental system of the regressive 
equation 0T71). Let f G Crd- Then the solution of the initial value problem 



Ly = f{t), y^'(to)=0, 0<2<n-l, 
is given hy y{t) = J^^y{t, s)f{s)As, where y{t, s) is the Cauchy function for 

It is known that y(t, s) := Hn-i(t, cr{s)) is the Cauchy function for y^" = 0, where Hn-i 
is a time scale generalized polynomial |4H Example 5.115]. The generalized polynomials 
Hk are the functions : — )> M, A; G Nq, defined recursively as follows: 

Ho{t,s)^l, H,+,{t,s)= [ Hk{r,s)Ar, /c = l,2,... 

J s 

for all s,t E T. If we let L[^{t, s) denote, for each fixed s, the derivative of Hk{t, s) with 
respect to t, then (cf. |1T| p. 38]) 

H^{t, s) = Hk-i{t, s) for A; G N, t G T^ 

Let a G M and h > 0, {h'Ij)a = {a, a + h,a + 2h, . . .}, and h = a + kh ioi some /c G M. 

From now on we restrict ourselves to the time scale T = {}iL)a, h > 0. Our main 
goal is to propose and develop a discrete-time fractional variational theory in T = (/iZ)^. 
We borrow the notations from the recent calculus of variations on time scales ISZlESlEn]- 
How to generalize our results to an arbitrary time scale T, with the graininess function fi 
depending on time, is not clear and remains a challenging question. 

We have a(t) = t+h, p(t) = t—h, fi(t) = h, and we will frequently write /"^(t) = f{a(t)). 
We put T = [a,b] n {hZ)a, so that = [a,p{b)] n {hZ)a and T''' = [a,p^{b)] n {hZ)a. 
The delta derivative coincides in this case with the forward /i-difference: 

nt) - fit) 



Pit) 
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li h = 1, then we have the usual discrete forward difference Af{t). 

The delta integral gives the /i-sum (or /i-integral) of /: / f(t)At = ''^^ f(kh)h. If we 

•J a ,_a 



have a function / of two variables, f(t, s), its partial forward /i-differences will be denoted 



by At^h and A^^^, respectively. We make use of the standard conventions Ylt=c /(^) ~ 0' 
c G Z, and Hi^^o /(^) ~ Oft^n, left fractional delta integration (resp., right fractional 
delta integration) of order u > is denoted by aA^'^f{t) (resp. tA^'"/(t)). Here, similarly 
as in Ross et al. |107] . where the authors omit the subscript t on the operator (the operator 
itself cannot depend on t), we write aAj^'^ f{t) (resp. ftA^'^/(t)). 

Before giving an explicit formula for the generalized polynomials on KL we introduce 
the following definition: 

Definition 71. For arbitrary x,?/ G M the /i-factorial function is defined by 



" ' r(f + l-y)' 
where F is the well-known Euler gamma function, and we use the convention that division 
at a pole yields zero. 

Remark 72. Before proposing Definition [7T] we have tried other possibilities. One that 
seemed most obvious is just to replace 1 by /i in formula (13. 2 p because if we do /i = 1 
on the time scale {h'L)a we have the time scale Z. This possibility did not reveal a good 
choice. 

Remark 73. For h = 1, and in accordance with the previous similar definition on Chapter [3l 
we write x^^-* to denote x^/^K 

Proposition 74. For the time-scale T = {h'L)a one has 

(t-s)^''^ 

Hk{t, s) := ^ — for all s,t e T and k e No . (4.4) 

K ! 

To prove (14. 4p we use the following technical lemma. Throughout this chapter the basic 
property (11.60 of the gamma function will be frequently used. 



Lemma 75. Let s G T. Then, for all t G T'' one has 

At,h 



(t-.)rn _(t-#) 



h 

{k + l)\ I k\ 
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Proof. The equality follows by direct computations: 



(t-.)rn _i rKt)-s)r) (t-.)r' 



{k + iy. j h[ {k + l)\ {k + l)\ 

h^+^ r T{{t + h- s)/h + l) T{{t-s)/h + l) 



h{k + 1)! i r((t + h- s)/h + l- (k + l)) T{{t - s)/h + 1 - (A; + 1)) 
/i^- ( iit-s)/h+l)T{{t-s)/h + l) T{{t-s)/h + l) 



{k + 1)! i {{t - s)/h - k)T{{t - s)/h - k) r((t - s)lh - k) 
\ T{{t-s)/h + l) 



k\ \T{{t- s)/h + l-k)] k\ 



Proof, (of Proposition 17^ We proceed by mathematical induction. For /c = 

- 0!^ r(^^ + i-o) - r(^^ + i) - ^ • 

Assume that fl4.4p holds for k replaced by m. Then by Lemma [731 



□ 



r — s 



H..»{t. s) = HUr, s)Ar = '-^^r , 

which is (I4.4p with k replaced by m + 1 . □ 

Let yi{t), . . . , Unit) be n linearly independent solutions of the linear homogeneous dy- 
namic equation y^" = 0. From Theorem [70] we know that the solution of (14. 3p (with 
L = A"- and to = a) is 

y{t)=A-f{t)= / ^ Z}"'!' fis)As = —- Y,{t-aikh)t-'^f{kh)h. 

^ ^ k=a/h 

Since y^'{a) = 0, z = 0, . . . , n — 1, then we can write that 

t/h—n 

^=<'/h (4.5) 

Note that function t — )■ (A~"/)(t) is defined for t = a + nh mod(/i) while function t — )■ /(t) 
is defined for t = a mod(/i). Extending (14.50 to any positive real value u, and having as an 
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analogy the continuous left and right fractional derivatives jST] , we define the left fractional 
/i-sum and the right fractional /i-sum as follows. We denote by J-r the set of all real valued 
functions defined on a given time scale T. 

Definition 76. Let a E R, h > 0, b = a + kh with k e N, and put T = [a,b] n {hZ)a. 
Consider / G Tj. The left and right fractional /i-sum of order u > are, respectively, the 
operators a^h" '■ J^t — ^ J^f+ and h^b" '■ -^t ~^ -^f^' = {t ± z/Zi : t G T}, defined by 

--1/ 

a^H^m = ^ - cris))t'^fis)As = Y^^Y.}^- aikh))t'^fikh)h 

,A^'^fit) = ^J^^Js-ait)Yr'^fis)As = ^^ J2 {kh-ait))t'^fikh)h. 



Remark 77. In Definition [TSl we are using summations with limits that are reals. For 
example, the summation that appears in the definition of operator aAy^" has the following 
meaning: 

J2 G{k) = G{a/h) + G{a/h + 1) + G{a/h + 2) + ■ ■ ■ + G{t/h - u), 
where t G {a + vh^ a + h + uh, a + 2h + uh, . . . ,a + kh +vh} with G M. 

6 

Lemma 78. Let u > be an arbitrary positive real number. For any t E T we have: (i) 
Xim^^^aAl^fit + vh) = fit); (ii) ^Al" f {t - vh) = fit). 

Proof. Since 



A-/(t + ^h) = —- I it + vh- ai.s))t'^fi.s)As 



l_Y,it + ^h-aikh))t~'^fikh)h 



h 

V 



it follows that Imiy^Q aA^^'^ f it + vh) = fit). The proof of (ii) is similar. □ 
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For any t eT and for any > we define aA°/(t) := /iA°/(t) := f{t) and write 



/J 

.A^/W = h-'fit - uh) + — — - / {s + ,^h- a{t))t'^f{s)As 

J- 1^^ + J-j Ja(t) 



Theorem 79. Let f E J-'j and > 0. For all t E T'^ we have 



To prove Theorem [72] we make use of a technical lemma: 



h 



r 



h 



r + 1) 



r ( + 1 - (z/ - 1) 



r( 



t+uh—s 
h 



r(^^ + i) r(^^ + 2) 

•p^ t+uh—s—h I j^^ 



+ (t+z./.-a(s));r^)/^(.) 



(4.6) 



(4.7) 



Lemma 80. Lei t G T*^. T/ie following equality holds for all s eT'^: 

{{t + uh-s)t'^f{s))) 

= {t + uh- a{s))t'^f%) -{v-l){t + uh- a{s))t'^f{s) . (4.8) 
Proof. Direct calculations give the intended result: 
A.,. ^(t + uh-s)t'^f{s)) 

= A,,, ((t + uh- s)t'^) fis) + it + uh- a{s))t'^ f^is) 



r(^^^±f^ + 1 - (z/ - 2)) 



{-{u-l)) + {t + uh-a{s))t'^f^{s) 



h 

= -{u -l){t + uh- a{s))t'^f{s) + {t + uh- a{s))\:~'^ f^s) 
where the first equality follows directly from (12.31) . 



□ 



Remark 81. Given an arbitrary t G it is easy to prove, in a similar way as in the proof 
of Lemma [HOl the following equality analogous to (14. 8 p : for all s G T'^ 



A.,. {{s + i.h-ait))t'^fis)i 

= {u-l){s + uh- a{t))t'^f-{s) + {s + uh- a{t))t'^fHs) • (4.9) 
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Proof, (of Theorem 179]) From Lemma |80] we obtain that 



{t + uh-s)t"'^f{s) 



s=t 



r(z/ + i)7a 



(4.10) 



J a 



We now show that {aAl''f{t + vh))^ equals (gJO]): 



{aAT:f{t + uh)Y 



h 



fait) 
J a 



V 



r(i. + i) 



/ir(z/ + 1) 



J a 

[a{t) + vh-a{s))t~^^f{s)As 



- r(t + i^/^-a(s))i'^-^V(=^)AJ +h'^~^ 

J a 



V 



r(^ + i)7a 



r Ai,, ((t + - a(s))i'^"')) ]{s)As + /i''-V/(t) 

J a 



□ 



Follows the counterpart of Theorem [79] for the right fractional h-sum: 
Theorem 82. Let f E J-'j and u > 0. For all t eT'^ we have 



KA-U\t-uh) 



V 



Proof. From (14. 9 p we obtain from integration by parts (item 2 of Lemma [27|) that 



(4.11) 



u{b + uh-a{t))t '^ 



r(^ + i)y. 



f(b) + h^f^{t)-uh^~'fia{t)) 



{u-l)is + uh-a{t))i^-'^f'^{s)As. 



(4.12) 



it) 
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We show that {h\''f{t - vh))^ equals fl412|l 



V 



hViu+l] 



h^]\t) + — — - / A,, ((. ^vh- a{t))t'A f{s)As - uh^-'f{a{t)) 



h^f^{t) 



V 



V 



b 



it) 



□ 



Definition 83. Let < a < 1 and set 7 := 1 — a. The left fractional difference aA^/(t) 
and the right fractional difference hA'^fit) of order a of a function f E J-'j are defined as 

aA-J{t) ■.= {^A-y {t + jh))^ and ,A^f{t):=-{,Arf{t-ih))^ 



for all teT^. 



4.3 Main Results 



Our aim is to introduce the /i-fractional problem of the calculus of variations and to 
prove corresponding necessary optimality conditions. In order to obtain an Euler-Lagrange 
type equation (cf. Theorem 1881) we first prove a fractional formula of /i-summation by parts. 



4.3.1 Fractional /i- summation by parts 

A big challenge was to discover a fractional /i-summation by parts formula within the 
time scale setting. Indeed, there is no clue of what such a formula should be. We found it 
eventually, making use of the following lemma. 
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Lemma 84. Let f and k be two functions defined on T'^ and T" , respectively, and g a 

2 

function defined on T'^ x T*^ . The following equality holds: 

fit) / g{t,s)k{s)As At= / k{t) I g{s,t)f{s)As At. 

Proof. Consider the matrices R = [/(a + h), f{a + 2h), ■ ■ ■ , f{b — h)], 

g{a + h, a)k{a) 

I g{a + 2h, a)k{a) + g{a + 2h, a + h)k{a + h) 



gib — h, a)k{a) 


+ g{b-h, 


a + h)k{a + h) + ■ 


■■ + g{b- 


h,b- 


2h)k{b 


-2h) 


g{a + h, a) 
















g{a + 2h, a) 




g{a + 2h, a + h) 


, Ck = 









gib - h,a) 




gib — h,a + h) 




9{b 


-h,b- 


2h) 



Direct calculations show that 



fit) 



g{t,s)k{s)As 



b/h~l 



At = ej2 f('^) E 9{^h,jh)k{jh) = h^R ■ Ci 

i=a/h j=a/h 



h^R ■ [k{a)C2 + k{a + h)C3 + ■ ■ ■ + k{b - 2h)Ck 



b/h-l b/h~l 

Ho-) ^ 9ijh,a)f{jh) + k{a + h) ^ g{jh,a + h)f{jh) 

j=a/h+l j=a/h+2 



b/h~l 



+ ... + k{b-2h) J2 9{jh,b-2h)f{jh) 

j=b/h-i 

b/h-2 b/h-1 (;,) , 

Y,H^h)h 9ijh,th)fijh)h= kit) / gi.s,t)fi.s)A.t 



i=a/h 



j=a{ih)/h 



At. 



□ 



Theorem 85 (fractional /i-summation by parts). Let f and g be real valued functions 
defined on T'^ and T, respectively. Fix < a < 1 and put 7 := 1 — a. Then, 

rp{b) 

f{t),Alg{t)At = hy{p{b))g{b) - hy{a)g{a) + / ,A°(,)/(t)<7'^(t)At 

J a 

^ :9ia) ( [\t + ^h- a)^r'^f{t)At - f {t + ^h- a(a))^^VW At) . (4.13) 



+ 



r(7 + i)- 
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Proof. By (14. 7p we can write 

J a J a 

= J' fit) [.A-^^^(t + ^h) + f^:^^it + lh- a)^r'^g{a) 

= r/(t)aA,"V(t+7^)At+ r 

J a J a 

Using (14. 6 p we get 

ff{t),A-''g^{t + jh)At 

J a 



At 



(4.14) 



it + ^h-a)^r'''f{t)gia)At. 



r(7 + i) 



m 



h^g^it) 



"Pib) 



At 



W / f{t)g^{t)At + 



7 



r(7 + i) 



/ g'^it) / {s + ^h-a{t))^r'''fis)AsAt 



Pib) 



= hy{pmg{b) - g{pm + / ^^(t),A-y (t - ^h)At, 

J a 

where the third equahty follows by Lemma [84l We proceed to develop the right hand side 
of the last equality as follows: 

Pib) 



rpyo) 

hV{p{h))[g{h) - g{p{h))] + / ^^(t),A-y(t - ^h)At 

J a 



t=p{b) 
t=a 



hV{p{h))[g{h) - g{pm + \g{t\A-l^f{t - ^h) 

rPib) 

g''mHA;;i)fit-^h))^At 

hV\p{b))g{b) - hy{a)g{a) 

o/ fPib) 

^ ^Ma) / (s + ^h- a{a))^r'^f{s)As + / (.A^"(,)/(i)) g'^{t)At, 



r(7 



where the first equality follows from Lemma [271 Putting this into fl4.14p we get fl4.13p . □ 

4.3.2 Necessary optimality conditions 

We begin to fix two arbitrary real numbers a and P such that a, /3 G (0, 1]. Further, 
we put 7 := 1 — a and z/ := 1 — /3. 
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Let a function L{t, u,v,w) : T'' x M x M x M — )■ M be given. We consider the problem of 
minimizing (or maximizing) a functional C : J^j — M subject to given boundary conditions: 

£(!/(■))= f L{t,y^{t),aKy{t).h^ly{t))^t^mm, y{a) = A, y{h) = B . (4.15) 

J a 

Our main aim is to derive necessary optimality conditions for problem (I4.15p . 

Definition 86. For / G J-r we define the norm 

11/11 = max ir (t)| + max UA°/(t)| + max |,Af /(t)|. 

A function y G J^f with y{a) = A and y{b) = B is called a local minimizer for problem 
f l4.15p provided there exists 6 > such that C{y) < C{y) for all y E J-'j with y{a) = A and 
y{b) = B and \\y — y\\ < 6. 

Definition 87. A function rj G J-j is called an admissible variation provided rj ^ and 
77(a) = 77(6) = 0. 

From now on we assume that the second-order partial derivatives L^^, L^^, L^y^, Lyy^, 
Lyy, and L„„„ exist and are continuous. 



First order optimality condition 

Next theorem gives a first order necessary condition for problem (I4.15p . i.e., an Euler- 
Lagrange type equation for the fractional /i-difference setting. 

Theorem 88 (The /;,- fractional Euler-Lagrange equation for problem ( I4.15p ). //y G J-'j is 
a local minimizer for problem f l4.15p . then the equality 

Ly[y]it) + /.A^(,)L4y](t) + aA^L^m = (4.16) 

holds for allt G T"^ with operator [■] defined by [y]{s) = (s, ?/'^(s), aA°y(s), sA^?/(s)). 

Proof. Suppose that y{-) is a local minimizer of £[■]. Let ?7(-) be an arbitrarily fixed 
admissible variation and define a function $ : (^~f^^\j IMon) ~^ ^ 

^s) = C[y{-) + evi-)]. (4.17) 

This function has a minimum at e = 0, so we must have $'(0) = 0, i.e., 

b 

L^iymv^it) + Ly[y]it)aKvit) + L^[y]{t)hA^bVit)\ = 0, 
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which we may write equivalently as 



Pib) 



hL^[y]{t)ri''{t)\t=,(^,) 



Using Theorem 15^ and the fact that 77(a) = T]{b) = 0, we get 



+ / 



^Afr/(t)At = 0. (4.18) 



Pib) 



LMit)aKriit)/^t -- 
for the third term in fl4.18p . Using (14.111) it follows that 

J a 



(4.19) 



{h+uh-a{t))t~'\{h) 



At 



UTjih) 



{h + vh-a{t))t-'^L4y]{t)At. 

(4.20) 



We now use Lemma [H] to get 

L^[y]{t)hA~ib)^''{t-vh)At 



V{u + l) 



Ja{t) 



At 



h'L^[y]{t)ri''{t)At 

"Pib) 



+ 



r(i^ + i) Ja 



/ L4y]{t) / {s + uh-ait))i'-'\^{s)As 



(4.21) 



At 



h''L4y]{t)v^it)At 



+ 



n^+l) J a 



rj^it) I {t + uh-a{s))i:-'^L^[y]{s)As 



At 



7i''{t)aAl^ {L^[y]){t + vh)At. 
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We apply again the time scale integration by parts formula ( Lemma l27|) . this time to f)4.21[l . 



to obtain, 

rb 

r/^(t),A-^(L^[y])(t + z//i)At 
rp{b) 

= / v'^it),A-''{L4y]){t + uh)At 



+ (T^ib) - vipmaA^,^ {L^[y]) {t + z//i)|i=,(,) 

J a 

+ r/(6)aA^'^ (t + i//i)|t=p(b) - r]{p{b))aA^'' {L^[y]) {t + i//i)|t=p(,) 

= r/(6),A^^ (t + i^/i)|t=p(fe) - r/(a)aA^'^ {L^[y]) [t + 

J a 

Since 77(a) = 77(6) = we obtain, from fl4.2ip and f l4.22p . that 

/ L^[y](t),A-^,)r/^(t)At = - / r/'^(t).Af (t)At , 

J a J a 

and after inserting in fl4.20p . that 

rb rp{b) 

/ L^[y]{t)hAlv{tW= / r/'^WaAf (L^[y])(t)At. (4.23) 

J a J a 

By (HTTg]) and <KTS^ we may write (H:TH]) as 

rp{b) 

/ L„[y](t) + ,A^(,) (L4y]) (t) + ,Af {L^[y]) (t) r/'^(t)At = . 

J a 

Since the values of rj^it) are arbitrary for t G T*^^, the Euler- Lagrange equation fl4.16p holds 
along y. □ 



The next result is a direct corollary of Theorem | 

Corollary 89 (The /i-Euler-Lagrange equation - cf., e.g., |37[[52]). Let T he the time scale 
KL, h > 0, with the forward jump operator a and the delta derivative A. Assume a,b eT, 
a <b. If y is a solution to the problem 

C{y{-))= I L(t,2/-(t),y^(t))At-^min, y{a) = A, y{b) = B , 

J a 

then the equality Lu{t,y''{t),y^{t)) - {L^{t,y''{t),y^{t)))^ = holds for all t E T''\ 
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Proof. Choose a = 1 and a L that does not depend on w in Theorem [881 D 
Remark 90. If we take /i = 1 in Corollary 189) we have that 

L,(t, y'^(t), Ay{t)) - ^L,{t, r{t), Ay(t)) = 

2 

holds for all t G T'' . This equation is usually called the discrete Euler- Lagrange equation, 
and can be found, e.g., in [TTj Chap. 8]. 



Natural boundary conditions 

If the initial condition y{a) = A is not present in problem fl4.15p (i.e., y{a) is free), 
besides the /i-fractional Euler-Lagrange equation (I4.16P the following supplementary con- 
dition must be fulfilled: 



-h^LJy]ia)+ ^ 



r(7 + i) 



[\t + ^h-aY^~'^LSmAt 

J a 

b 

it + ^h-a{a))\;'-'^L4y]{t)^t] + L4y]{a) = 0. (4.24) 

a(a) 



Similarly, if y{b) = B is not present in (I4.15P {y{b) is free), the extra condition 
hL^Mb)) + h^LMipib)) - h^L^,[y]{pib)) 



V 

+ 



- I'^Pi^) + '^h- a{t))[:-'^L^m)^t^ = (4.25) 

is added to Theorem |88l We leave the proof of the natural boundary conditions fl4.24p and 
(I4.25P to the reader. We just note here that the first term in fl4.25p arises from the first 
term of the left hand side of (I4.18P . 



Second order optimality condition 

We now obtain a second order necessary condition for problem fl4.15p . i.e., we prove a 
Legendre optimality type condition for the fractional /i-difference setting. 
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Theorem 91 (The /i-fractional Legendre necessary condition). If y E J-x is a local mini- 
mizer for problem f l4.15|) . then the inequality 

h'^Luuim) + 2/i^+^L„4y](t) + 2h''+\u - l)L^4m) + ^'^(7 - lfL,,[y]{a{t)) 
+ 2/1^+^(7 - l)L,^[y\{a{t)) + 2h''+\v - 1)L, 



h^^u-l?L 



+ h'-L^M{a{t))+ / h'L. 

J a 

+ h^L,,[y]{t)+ [ h^L„ 

2 

holds for all t gT'^ , where 



uil — y) , , , , , 



WW 

2 



{u-2) 
h 



As 

2 



(4.26) 



^)(^|^(^ + 7/^-^Kt)))r'^) As>0 
(t,rW,aA,"y(t),tAf^(t)). 



Proof. By the hypothesis of the theorem, and letting $ be as in fl4.17p . we have as necessary 
optimahty condition that $"(0) > for an arbitrary admissible variation ?7(-). Inequality 
$"(0) > is equivalent to 



Luuiymiv^it)? + 2L^Mit)v''{t)aA'',,v{t) + 2L„^[y](t)r/-(t),Afr/(t) 
+L™[y](t)(,A^r/(t))2 + 2L,4y]{t)aA'^7]{t)HA^,7]{t) + L^^{t){HA^,r]{t))^] At > 0. (4.27) 



2 I /i if t = cr(r); 

Let r G T*^ be arbitrary, and choose 77 : T — > M given by 77 (t) = < 

[ otherwise. 

follows that 77(a) = 77(6) = 0, i.e., 7/ is an admissible variation. Using fl4.7p we get 



[Luulymiv^ityf + 2L^Mit)v''{t)aA'^HV{t) + L^SmiaA-^vit)?] At 

+ 2LUy]{tW{t) [h^V^it) + ^T^:^ J\t + ^h- a{s))t'Y{s)As^ 



At 



/7^L„4y](r) + 2h'<^^L^,[y]{r) + K<+^ L,,[y]{T) 



+ / Lyy 

lair) 



h-'v^it) + 



7 



r(7 + 



-j\t + ^h-ais))t'^vHs)As^ At. 



It 
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Observe that 



/.2^+i(7-l)%.[y]Kr)) 

rb 



1 



r(7 + i) 



j\t + ^h-ais))^^-'^7]^is)As^ At 



Let t G [a^{T),p{b)] n hZ. Since 



7 



r(7 + i) 



(t + ^/,_^(,))(7-V(^)A. 

Mr) 



7 



r(7 + i) 



7 



r(7 + i) L 



r(7 + i) 



Ja{r) 

(¥ + 7-i)r(^^ + 7-i)-(^)r(^^ + 7-i) 



(4.28) 



:¥)r(¥) 



! 7(7- 1) 

r(7 + i) 



we conclude that 



a2(r) 



7 



r(7 + i) 



j\t + ^h-a{s))^^~'^r]^{s)As^ At 

7(7 - 1) 



b 



Note that we can write tAl'r]{t) = —hA^^^^7]'^{t — uh) because 77(6) = 0. It is not difficuh 
to see that the following equality holds: 



2LuM{t)T^''{t)hAlT^{t)At = - 2Lu4y]{t)v''{t)hA-/,^7]^{t - uh)At 

J a 

= 2h^+''Lu4y]{T){u-l). 
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Moreover, 

rb 

2L,4y]{t)aKv{t)h^h{t)^t 
= -2 f L,^[y]{t){{h''r]\t) + 



1 



r(7 + l) J a 



(t + 7/i-a(s))^^-V(^)As 



h'ri^it) + 



V 



{s + uh-a{t)t-'^^\s)^s 



r(i^ + l)^(t) 

= 2/i^+-+i(z/ - l)L,„^[y](r) + 2/1^+^^+1(7 - l)L,^[y](a(r)). 
Finally, we have that 

L^^[y](t)(,Afr/(t))2At 



At 



r(i. + i) 



f {s + uh-a{t)ir'^r^^{s)^s 

Ja{t) 



At 



■/ (s + z//i-(T(t))i""'^r/^(s)As At 



(r + z./.-a(t))!r^)-(a(r) + z./.-a(t))r 



Similarly as we did in (I4.28p . we can prove that 



(^-1) 



r(^ + i) 



{T + uh- - (a(r) + uh - a 



z/(l - z/) 



(1.-2) 



Thus, we have that inequality (I4.27P is equivalent to 



h\h^L^S]{t)+2h^^'L, 



a(<T(t)) Vr(7 + 1) 



As 



+ 2h''+^L^^[y]{t){u - 1) + 2W+''{u - l)L,^[y]{t) 
+ 2/1^+^^(7 - + - vf + /i2'^L^,^[y](a(t)) 

V(l-z.) 



+ / h-'L, 



{t + uh - a{s)Y"' ] AsS>0. (4.29) 
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Because h > 0, fl4.29p is equivalent to fl4.26p . The theorem is proved. □ 
The next result is a simple corollary of Theorem [HH 

Corollary 92 (The /i-Legendre necessary condition - cf. Result 1.3 of [37]). Let T be 

the time scale KL, h > 0, with the forward jump operator a and the delta derivative A. 
Assume a,b & T , a < b. If y is a solution to the problem 

nb 

C{y{-))= / L(t,y'^(t),y^(t))At ^min, y{a) = A, y{b) = B , 

J a 

then the inequality 

h^L^S]{t) + 2hL^,[y]{t) + L,,[y]{t) + L,,[y]{(r{t)) > (4.30) 
holds for allte T"-', where [y]{t) = {t,y''{t),y^{t)) . 

Proof. Choose a = 1 and a Lagrangian L that does not depend on w. Then, 7 = and 
the result follows immediately from Theorem [9T] □ 

Remark 93. When h goes to zero we have a{t) = t and inequality (I4.30p coincides with 
Legendre's classical necessary optimality condition Lyy[y](t) > (cf., e.g., |114| ). 



4.4 Examples 

In this section we present some illustrative examples. 
Example 94. Let us consider the following problem: 

^(y) = ll (oA||/(t))'At-^min, y(0) = , y(l) = 1 . (4.31) 

We consider fl4.3ip with different values of h. Numerical results show that when h tends to 
zero the /i-fractional Euler-Lagrange extremal tends to the fractional continuous extremal: 

when h ^ (I4.3ip tends to the fractional continuous variational problem in the Riemann- 
Liouville sense studied in [3l Example 1], with solution given by 

yit) = I f r • (4.32) 

^ Jo [{l-x){t-x)Y 

This is illustrated in Figure 14.11 

In this example for each value of h there is a unique ^-fractional Euler-Lagrange extremal, 

solution of (I4.16p . which always verifies the /i-fractional Legendre necessary condition 
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a 



m 



1.0 




0.9 



0.8 



0.7 



0.6 



0.5 



0.4 



0.3 



0.2 



0.1 











Figure 4.1: Extremal y{t) for problem of Example [Ml with different values of h: h = 0.50 
(•); h = 0.125 (+); h = 0.0625 (*); /i = 1/30 (x). The continuous line represent function 



We begin by considering problem fl4.33p with a fixed value for a and different values of h. 
The extremals y are obtained using our Euler-Lagrange equation f l4.16p . As in Example [Ml 
the numerical results show that when h tends to zero the extremal of the problem tends to 
the extremal of the corresponding continuous fractional problem of the calculus of variations 
in the Riemann-Liouville sense. More precisely, when h approximates zero problem (I4.33P 
tends to the fractional continuous problem studied in [HI Example 2]. For a = 1 and h ^ 
the extremal of fl4.33p is given hj y(t) = ^t{l — t), which coincides with the extremal of 
the classical problem of the calculus of variations 



(1132]). 



Example 95. Let us consider the following problem: 




(4.33) 




This is illustrated in Figure 14.21 for h = i, i = 1,2, 3, 4. 
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m 

0.15 




Figure 4.2: Extremal y{t) for problem fl4.33p with a = 1 and different values oi h: h = 0.5 
(•); h = 0.25 (x); h = 0.125 (+); h = 0.0625 (*). 

In this example, for each value of a and h, we only have one extremal (we only have 
one solution to (14. 16p for each a and h). Our Legendre condition fl4.26p is always verified 
along the extremals. Figure 14.31 shows the extremals of problem fl4.33p for a fixed value of 
h [h = 1/20) and different values of a. The numerical results show that when a tends to 
one the extremal tends to the solution of the classical (integer order) discrete-time problem. 




Figure 4.3: Extremal y{t) for fl4.33p with h = 0.05 and different values of a: a = 0.70 (•); 
a = 0.75 (x); a = 0.95 (+); a = 0.99 (*). The continuous line is y{t) = it(l - 1). 
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Our last example shows that the /i-fractional Legendre necessary optimality condition 
can be a very useful tool. In Example [961 we consider a problem for which the /i-fractional 
Euler-Lagrange equation gives several candidates but just a few of them verify the Legendre 
condition ( K^ . 

Example 96. Let us consider the following problem: 



£{y)= UAty{t)f + e{,A^y{t)fAt^mm, y{a) = , y{b) = 1 . (4.34) 



For a = 0.8, (3 = 0.5, h = 0.25, a = 0, 6 = 1, and 9 = 1, problem fl4.34p has eight different 
Euler-Lagrange extremals. As we can see on Table 14.11 only two of the candidates verify 
the Legendre condition. To determine the best candidate we compare the values of the 
functional C along the two good candidates. The extremal we are looking for is given by 
the candidate number five on Table 14. 1[ 



# 






HI) 




Legendre condition (]4.26p 


1 


-0.5511786 


0.0515282 


0.5133134 


9.3035911 


Not verified 


2 


0.2669091 


0.4878808 


0.7151924 


2.0084203 


Verified 


3 


-2.6745703 


0.5599360 


-2.6730125 


698.4443232 


Not verified 


4 


0.5789976 


1.0701515 


0.1840377 


12.5174960 


Not verified 


5 


1.0306820 


1.8920322 


2.7429222 


-32.7189756 


Verified 


6 


0.5087946 


-0.1861431 


0.4489196 


10.6730959 


Not verified 


7 


4.0583690 


-1.0299054 


-5.0030989 


2451.7637948 


Not verified 


8 


-1.7436106 


-3.1898449 


-0.8850511 


238.6120299 


Not verified 



Table 4.1: There exist 8 Euler-Lagrange extremals for problem (I4.34p with a = 0.8, f3 = 0.5, 
h = 0.25, a = 0, b = 1, and 6 = 1, but only 2 of them satisfy the fractional Legendre 
condition (K26\i . 

For problem fl4.34p with a = 0.3, h = 0.1, a = 0, b = 0.5, and 6' = 0, we obtain the 
results of Table 14.21 there exist sixteen Euler-Lagrange extremals but only one satisfy the 
fractional Legendre condition. The extremal we are looking for is given by the candidate 
number six on Table 14.21 
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# 


y(o.i) 


y(o.2) 


y(0.3) 


y{OA) 




(|4.26|) 


1 


-0.305570704 


-0.428093486 


0.223708338 


0.480549114 


12.25396166 


No 


2 


-0.427934654 


-0.599520948 


0.313290997 


-0.661831134 


156.2317667 


No 


3 


0.284152257 


-0.227595659 


0.318847274 


0.531827387 


8.669645848 


No 


4 


-0.277642565 


0.222381632 


0.386666793 


0.555841555 


6.993518478 


No 


5 


0.387074742 


-0.310032839 


0.434336603 


-0.482903047 


110.7912605 


No 


6 


0.259846344 


0.364035314 


0.463222456 


0.597907505 


5.104389191 


Yes 


7 


-0.375094681 


0.300437245 


0.522386246 


-0.419053781 


93.95316858 


No 


8 


0.343327771 


0.480989769 


0.61204299 


-0.280908953 


69.23497954 


No 


9 


0.297792192 


0.417196073 


-0.218013689 


0.460556635 


14.12227593 


No 


10 


0.41283304 


0.578364133 


-0.302235104 


-0.649232892 


157.8272685 


No 


11 


-0.321401682 


0.257431098 


-0.360644857 


0.400971272 


19.87468886 


No 


12 


0.330157414 


-0.264444122 


-0.459803086 


0.368850105 


24.84475504 


No 


13 


-0.459640837 


0.368155651 


-0.515763025 


-0.860276767 


224.9964788 


No 


14 


-0.359429958 


-0.50354835 


-0.640748011 


0.294083676 


34.43515839 


No 


15 


0.477760586 


-0.382668914 


-0.66536683 


-0.956478654 


263.3075289 


No 


16 


-0.541587541 


-0.758744525 


-0.965476394 


-1.246195157 


392.9592508 


No 



Table 4.2: There exist 16 Euler-Lagrange extremals for problem (I4.34p with a = 0.3, 
h = 0.1, a = 0, b = 0.5, and ^ = 0, but only 1 (candidate ^6) satisfy the fractional 
Legendre condition fl4.26p . 
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4.5 Conclusion 

In this chapter we introduce a new fractional difference variational calculus in the 
time-scale (/iZ)^, h > and a a real number, for Lagrangians depending on left and right 
discrete-time fractional derivatives. Our objective was to introduce the concept of left and 
right fractional sum/difference (cf. Definition 1761) and to develop the theory of fractional 
difference calculus. An Euler-Lagrange type equation fl4.16p . fractional natural boundary 
conditions ( I4.24p and f l4.25p . and a second order Legendre type necessary optimality con- 
dition (I4.26p . were obtained. The results are based on a new discrete fractional summation 
by parts formula (I4.13P for {hZ)a- Obtained first and second order necessary optimality 
conditions were implemented computationally in the computer algebra systems Maple and 
Maxima (the Maxima code is found in Appendix iBl). Our numerical results show that so- 
lutions to the considered fractional problems become the classical discrete-time solutions 
when the fractional order of the discrete-derivatives are integer values, and that they con- 
verge to the fractional continuous-time solutions when h tends to zero. Our Legendre type 
condition is useful to eliminate false candidates identified via the Euler-Lagrange fractional 
equation. The results of the chapter are formulated using standard notations of the the- 
ory of time scales | ^ [69 | [76] because we keep in our mind the desire to generalize the 
present results to an arbitrary time scale T. Undoubtedly, much remains to be done in 
the development of the theory of discrete fractional calculus of variations in (/iZ)^ here 
initiated. 

4.6 State of the Art 

The results of this chapter are published in [3l] and were presented by the author at 
the FSS'09, Symposium on Fractional Signals and Systems, Lisbon, Portugal, November 
4-6, 2009, in a contributed talk entitled The fractional difference calculus of variations. 

Recently, Ferreira and Torres [53] continued the development of the theory presented 
here with handy tools for the explicit solution of discrete equations involving left and right 
fractional difference operators. 
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Chapter 5 



Fractional Derivatives and Integrals on 
arbitrary T 



In this chapter we introduce a fractional calculus on time scales using the theory of 
delta dynamic equations. The basic notions of fractional order integral and fractional order 
derivative on an arbitrary time scale are proposed, using the inverse Laplace transform on 
time scales. Useful properties of the new fractional operators are proved. 



Recently, two attempts have been made to provide a general definition of fractional 
derivative on an arbitrary time scale (see [T8 |I104] ). These two works address a very 
interesting question, but unfortunately there is a small inconsistency in the very beginning 
of both studies. Indeed, investigations [T8 |I104| are based on the following definition of 
generalized polynomials on time scales /iq, : T x T — M: 



Although the recurrence formula above address a very interesting question about how to 
define fractional generalized polynomials on time scales, recursion (15. ip provides a definition 
only in the case a G Nq, and there is no hope to define polynomials for real or complex 
indices a with (15.11) . 



5.1 Introduction 



ho{t,s) = 1, 




(5.1) 
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Here we propose a different approach to provide a general definition of fractional deriva- 
tive on an arbitrary time scale based on the Laplace transform |42] . 

The chapter is organized as follows. In Section 15.21 we review the basic notions of 
Laplace transform on M (Section 15.2. ip and we introduce some necessary tools from time 
scales (Section I5.2.2p . Our results are then given in Section [5751 we introduce the concept 
of fractional integral and fractional derivative on an arbitrary time scale T (Section 15. 3. ip ; 
we then prove some important properties of the fractional integrals and derivatives (Sec- 
tion [532]). 

5.2 Preliminaries 

We start this section recalling results on classical Laplace transforms. All necessary 
concepts to understand the Laplace transform definition on time scales are supplied. 

5.2.1 Laplace transform on M as motivation 

In order to motivate our idea we start with some concepts of continuous classical Laplace 
transform. 

Definition 97. The Laplace transform of / : M — )■ M is defined by 



where 2; G C is chosen so that the integral converges absolutely. 
Remark 98. It is well known that if 

1. / is piecewise continuous on the interval < t < A ioi any positive A; 

2. 1/(^)1 < Me"* when t >T, for any real constant a and any positive constants M and 
T (this means that / is of exponential order, i.e., its rate of growth is not faster than 
that of exponential functions); 

then the Laplace transform (15. 2p exists for z > a. 




(5.2) 
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Because of the usefulness of the Laplace transform of derivatives in our work, we now 
refer what happens in the continuous case, using the following theorem that can be easily 
proved using the definition of Laplace transform and integration by parts. 

Theorem 99. Suppose f is of exponential order, and that f is continuous and f is piecewise 
continuous on any interval < t < A. Then 

C[f'] = zC[f]{z)-f{0) . (5.3) 

Remark 100. Applying the theorem multiple times yield 

= 2-£[/]_2-i/(o)-^-2/'(0)-. . .-^V^""'n0)-2/^""'n0)-/^""'H0). (5.4) 

The above equality is an extremely useful tool of the Laplace transform for solving linear 
ODEs with constant coefficients because it converts linear differential equations to linear 
algebraic equations that can be solved easily. 

Next proposition gives the Laplace transform of the Caputo fractional derivative and 
was an inspiration for our work. 

Proposition 101. f7^ Let a > 0, n be the integer such that n — 1 < a < n, and 

f a function satisfying f G C"(M+), G Li{0,ti), ti > 0, and |/(")(^)l < Bei°\ 

t > ti > 0. If the Laplace transforms C[f]{z) and C[f^"'']{z) exist, and lim f^^\t) = for 
k = 0, . . . ,n — 1, then 

n-l 

^ [^D:f] {z) = z-c [f] {z) - J2 f^'\^)z^-'-' . 

k=0 

Remark 102. If a G (0, 1], then C [^D^f] (z) = z'^C [/] (z) - f{0)z"-\ 

5.2.2 The Laplace transform on time scales 

Looking to Definition [27] we can observe that, in classical calculus, the exponential 
function assumes an important role in the theory of Laplace transforms. 

Following the same direction, we would like to have a function on time scales that serve 
the same purpose as the exponential function does in the real case. 

Before we can define such function, we need several concepts in order for the definition 
to make sense. The first concept that we need is the concept of regressivity. 
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Definition 103. A function p : T'' — )■ M is regressive provided 

holds for all t G T"^. We denote by TZ the set of all regressive and rd-continuous functions. 
The set of all positively regressive functions is defined by 

7^+ = {p e 7^ : 1 + fi{t)p{t) > 0, for all t E T^}. 

Definition 104. The function (ep)(t) = - 1+ J(t)p(f) for all teT^'andpe U. 

Now we define the exponential function on time scales (also called, in some literature, 
the generalized exponential function) as the solution of the IVP (15. 5p . 

Theorem 105. |^ Theorem 1.37] Suppose p eTZ and fix to G T. Then the initial value 
problem 

= Pit)y, yito) = 1 (5.5) 

has a unique solution on T. 

The solution of the IVP (15. 5p is the exponential function on time scales, and is given 
by ep{-,to), where 

epit,s) = exp{ ^,,^rMr))Ar\ with U^) = { , ^ • (5-6) 

Us J [ z it h = 

In (15. 6p C,h is the cylinder transformation defined for /i > from the set : = 
{z eC: -i} to the set Z?, := {2 G C : -f < Im{z) < f }. The set Ch is called the 
Hilger complex plane |41) . 

Remark 106. For h = 0, Cq := C. 

Definition 107. [HI p. 52] Let h > and z E Ch- The Hilger real part of z is defined by 

, ^ \zh + 1\ — 1 , ^ 

Rchiz) := J (5.7) 

and the Hilger imaginary part of z by 

, , Arqizh + 1) , , 

Imh{z) := ^, (5.8) 

where Arg{z) denotes, as usual, the principal argument of z. 
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Remark 108. If we consider a time scale with graininess function fi{t) > and a z belonging 
to := G C : z 7^ replace h by in the right-hand sides of (15. 7p 

and (15. 8p and represent them, respectively, by Re^{z) and Im^{z). 

Now we recall some properties of the exponential function stated in Theorem 2.36 

of nu. 

Theorem 109. giy IfpeTZ, then, 

1. eQ(t,s) = 1 and ep{t,t) = 1; 

2. ep{a{t), s) = (1 + fi{t)p{t))ep{t, s); 

^- t;^^ = eQp{t, s) . 

Remark 110. Item 2 of Theorem 11091 has special importance because it provides the ability 
to rewrite ep{a{t),s) without the symbol cr(t) and on that way simplify expressions of 
exponential function of elementary functions and by consequence the Laplace transform. 

With all previous necessary concepts already given we are in conditions to understand 
the Laplace transform on time scales and provide some of their properties which are needed 
across this chapter. 

The main results presented in this section were taken from |32] and/or |47) . 

Definition 111. j32] We define the generalized Laplace transform of a regulated function 
/ : T — )■ C, where T denotes a time scale which is unbounded above and contains zero, by 

POO 

Cj[f]{z) = F{z) := f{t)eeMt),0)^t for z e D[fl (5.9) 
Jo 

where D[f] consists of all 2; G C for which the improper integral exists and for which 
1 + ^{t)z ^ for all t G T. 

Remark 112. In view of Definition II IH the Laplace transform C of Proposition 1101 1 can be 
written as C^. 

Throughout the rest of chapter, T is an arbitrary time scale with bounded graininess, 
i.e., < fimin < /^(^) < /^maz ioT all t G T. Let to ^ T be fixed. 

Definition 113. [1^ The function / : T — )■ R is said to be of exponential type I if 
there exist constants M, c > such that |/(t)| < Me'^*. Furthermore, / is said to be of 
exponential type II if there exist constants M, c > such that \ f{t)\ < Medt, 0). 
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The time scale exponential function ec(t, 0) is of type II while generalized polynomials 
hk{t, 0) are of type I. 

Theorem 114. [Jl^ If f is of exponential type II with exponential constant c, then the 
delta integral f(t)eQz{o'{t),0)At converges absolutely for z G D. 

Similarly to the classical Laplace transform (cf. formula (15. 3p ) we can state the following 
theorem: 

Theorem 115. fj^ Theorem 1.2] If f : T ^ C is such that f^ is regulated, then 

Cj[f'']{z) = zCj[f]{z)-f{0) (5.10) 
for all z G T>[f] such that limf_>,oo {f{t)eQzit)) = 0. 

Remark 116. To prove Theorem 11151 we first make use of the Laplace transform, next 
apply integration by parts formula, and then use some of the properties of the exponential 
function mentioned before to rewrite the result and get what the theorem claims. 

Using the last result and mathematical induction one obtains the following result: 

Proposition 117. Let / : T — > R 6e such that is regulated and F represents the 
generalized Laplace transform of f . Then, 

n-l 

CtU^z) = z-F{z) - ^^"-'=-V^'^(0) (5.11) 

fc=0 

for all z G V[f] such that limf^oo (/'^'(^)eez(^)) = (0 <k <n). 

Proof. The proof will be done by mathematical induction on n. Before starting the induc- 
tion process it's important to note that the fact that f^" is regulated implies that all delta 
derivatives with order below n are also regulated. We know from Theorem 11151 that 

Cr[f^]{z) = zCr[f]iz)-f{0). 
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Assuming that (15. lip is true we have 

£t[/^""](^) = zCrlf'liz) - no) 

= z (z-F{z) - ^^^-'-v^^o)) - no) 

\ k=o J 

n-l 

= z-^'F{z)-J2z-~'f'^\0)-z'f^\0) 

k=0 
n 

fc=0 

□ 

Using Theorem I118[ Davis et al. |17] estabhshed Theorem 11191 below. 

Theorem 118. fj^ Theorem 1.2] Let F denote the generalized transform for / : T — )■ M. 
Then: 

L F{z) is analytic in Re^{z) > Re^{c); 

2. F{z) is hounded in Re^{z) > Re^{c); 

3. hm|,|^ooi"(^) = 0. 

Since our idea is to apply the Laplace transform theory on time scales, particularly the 
inverse of Laplace transform on time scales, to define the fractional integral (cf. Defini- 
tion [T2S]) and fractional derivative (cf. Definition I126p on time scales. Theorem 11191 plays 
an important role in our work by providing sufficient conditions for the existence of such 
inverse and a formula to evaluate it. 

Theorem 119 (Inversion formula of the Laplace transform |47|). Suppose that F is analytic 
in the region Re^{z) > Re^{c) and F{z) — )■ uniformly as \z\ — > oo in this region. 
Assume F has finitely many regressive poles of finite order {zi, 2:2, ... , Zn} and F-^{z) is 
the transform of the function f(t) on M that corresponds to the transform F{z) = Fj{z) of 
fit) on T. If 



pc+ioo 








Mz) 


\dz\ 


1 c—ioo 







then 

n 

f{t) = Y,R^s,=zMt,o)F{z) 

i=l 

has transform F{z) for all z with Re{z) > c. 
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Example 120. [HI Example 1.1] Suppose F{z) = ^. Then 

e.(t,0) 



C^^[F]{z)=f{t)=Res,=o- 



^2 



Applying Theorem 11191 for F[z) = ^ we get 

C:,^[F]{z) = f{t) = h,{t,Q). 

Remark 121. Using induction arguments over the results of Example 11201 it's easy to prove 
that the inverse transform of F[z) = jw+t, G N, is hk{-, 0). 

Remark 122. The inversion formula also gives the claimed inverses for any of the elementary 
functions that were presented in the table of Laplace transform in |42| . 

The following lemma establish a relation between the solution of (15. 5p with power series 
on time scales. 

Lemma 123. f39i Lemma 4-41 ^^r all z E C and t E T with t > a, the initial value 
problem 

y^ = zy, y{a) = l (5.12) 
has a unique solution y that is represented in the form 

+ 00 

y{t) = Y,z''hk{t,a) (5.13) 

k=0 

and satisfies the inequality 

1^(^)1 < el-K*-"). (5.14) 

Remark 124. Let a = 0. If we use jointly Lemma 11231 and the well-known fact that the 
unique solution of IVP ( 15. 5 p is ez{t, 0) we have the equality 



e,{t,0) = J2^'hk{t,0) 



k=0 



that allows us to see that generalized polynomials hk{t,0) is the reason to have different 
inverse Laplace images when we have different time scales. 



78 



5.3. MAIN RESULTS 



5.3 Main Results 

We begin by introducing the definition of fractional integral and fractional derivative 
on an arbitrary time scale T. 

5.3.1 Fractional derivative and integral on time scales 

Simultaneously, the generalized Laplace transform on time scales gives unification and 
extension of the classical results. Important to us, the Laplace transform of the A-derivative 
is given by the formula Cj[f^]{z) = zF{z) — /(O). Our idea is to define the fractional 
derivative on time scales via the inverse Laplace transform formula for the complex function 
G{z) = z"' Cj[f]{z) — /(0"*")2;°~^. Furthermore, for a G {n — l,n], n G N, we use a 
generalization of ( 15.1 ip to define fractional derivatives on times scales for higher orders a. 

Definition 125 (Fractional integral on time scales). Let a > 0, T be a time scale, and 
/ : T — M. The fractional integral of / of order a on the time scale T, denoted by Ijf, is 
defined by 



F(z) 



it) 



z" 

Definition 126 (Fractional derivative on time scales). Let T be a time scale, F{z) = 
Cf[f]{z), and a G (n — l,n], n G N. The fractional derivative of function / of order a on 
the time scale T, denoted by /"-"^ is defined by 

n-l 



(t). (5.15) 



'F(z) -^/^'(0+)z"-'=-^ 

fc=0 

Remark 127. For a G (0, 1] we have 

/(°)(t) = £^^ [z-F{z)-fiO^)z--'] (t). 
Moreover, if we use T = M we have 

/(")(t) = C^' [z'^Fiz) - /(0+)."-i] (t) = C^' [C[^,D2f]{z)] (t) = {^D:f) (t) . 

5.3.2 Properties 

We begin with two trivial but important remarks about the fractional integral and the 
fractional derivative operators just introduced. 
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Remark 128. As the inverse Laplace transform is linear, we also have linearity for the new 
fractional integral and derivative: 

I!^iaf + bg){t)=ai;^fit) + bi;^git), 
[af + hgf'^\t) = af^\t) + hg^'^\t). 

Remark 129. Since /lo(^) = 1 for any time scale T, from the definition of Laplace transform 
and the fractional derivative we conclude that hl^\t) = 0. For a > 1 one has h[''\t) = 0. 

We now prove several important properties of the fractional integrals and fractional 
derivatives on arbitrary time scales. 

Proposition 130. Let a E {n — l,n], n E N. If k <n — I, then 
Proof. From fl5.15p it follows that 



hi"\t,0) = C', 



yk+l 



(t) = o. 



i=0 



it) = 



yk + l 



J2hk^^iO)z''-'-' 



□ 



Proposition 131. Let a E {n — l,n\, n E N. If k > n, then 



1 



yk+l-a 



it) 



Proof. From (I5.15P we have 



h^k\t)=C^' 



n-l 



,a—i—l 



i=0 



yk + l 



it) = C:^' 



z 



k+l-a 



it) = 
it). 



n-l 



yk+l 



1=0 



it) 



□ 



Proposition 132. Let a E (n — l,n], n E N. If c{t) = m, m E M, then 
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Proof. From the linearity of the inverse Laplace transform and the fact that 
it follows that 

c(")(t) = (m- = (m- /io(t,0))(") = m- h''^\t,0) = m-0 = 0. 

Proposition 133. Let a,(3 > 0. Then, 

(T^f) (t) = I-^'fit) . 

Proof. By definition we have 

4 (I^f) it) = C^' [z~^C^ [/^/]] it) = C^' [s~-~^F{z)] it) = /^^/(t) . 

Proposition 134. Leta,f3 G (0,1]. 

(a) Ifa + (3< 1, then 

(b) Ifl<a + (3<2, then 

^fia^) ^ ^(.+,)^^) _ ^-1 [/-i/W(o)] (t) + [z'^^^-'f^iO)] (t) 

(c) If 13 ^ (0,1], then 

Proof, (a) Let a + /3 < 1. Then 

(/("))^^^ (t) = [.^/: (.) - .^-V^^Ho)] (t) 

= [z"+/^F(z) - ^°+^-V(0)] it) - C^' [^^^V^-HO)] (t) 

= f(»^p\t) - [2^-V^"Ho)] (t). 
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(b) For l<a + /3<2we have 

fi»+P)^t) = C^^ [z'^+fC [f] (z) - z"+^~V(0) - ^"+''-^'^(0)] (t) . 

Hence, 

i&Y^ (t) = f^'^^'Kt) + c^' [z-+'-'f^{o)] (t) - c^' [z^-'f^-Ho)] (t) . 

(c) The intended relation follows from (b) by choosing a = 1. □ 
In general (/'•"'') '''^^ (t) (Z^^-')*""'' (t)- However, the equality holds in particular cases: 

Proposition 135. If a + /S < 1 and /(O) = 0, then ^'^^ (t) = (/(^^)^°^ (t). 

Proof. It follows from item (a) of Proposition 11341 □ 

Proposition 136. Let a G {n — l,n], n G N, and lim f^'^it) = f^'^{0^) exist, k = 
0, . . . ,n — 1. The following equality holds: 

n-l 

i^ if^-^) it) = fit) -j2f^\o^)h,it) . 

Proof Let Fiz) = Cj [/] iz). Then, 



fc=0 



I^ (/(")) it) = C^' [z-Cr (.)] it) = 



^a-k-l 

\k , ^ Z 



Fiz)-J2f^\^- 



n-l 



k=0 



k=0 
n-l 



it) 



fit)-Y,f^\0^)^T' 3TT (t) = /(t)-5^/^''(0+)/.,(t). 



k=0 



□ 



Proposition 137. Let a e in - l,n], n e N, and Fiz) = Cj [f] iz). If lim Fiz) = and 

z—^oo 

lim 1^ = 0, A; G {1, ... , n}, then (/^/)(") (t) = fit). 



Proof. Firstly let us notice that lim iljf) it) = for k G {0, . . . , — 1}. For that 



we check that lim zC [Ijf] iz) = lim z^^ = 0. Nextly let us assume that it holds for 



i = A; — 1 and let us calculate 



lim zC 


'ilTf)^[ 


iz) = lim z 


I: 


z—^oo 


2— >00 








= lim z 


J. 
















82 



km 



k-l 



i=0 



2— i-OO Z' 



a—k—1 



5.4. STATE OF THE ART 



Then we easily conclude that 



it) = fit) . 



□ 



The convolution of two functions / : T — t- M and (7 : T x T — M on time scales, where 
g is rd-continuous with respect to the first variable, is defined in [IHIHT]: 

if*9)it)= f fiT)g{t,a{T))AT. 
Jo 

As function g we can consider, e.g., ec(t,to) or hkit^to). 
Proposition 138. LettQ e T. If a e (0,1), then 

if * 9i; to))^"^ it) = (/(^) * g^'^K; to)) W = (/^"^ * 9i; to)) it) , 
where we assume the existence of the involved derivatives. 

Proof. From the convolution theorem for the generalized Laplace transform Theorem 
2.1], 



(/*<7(-,to))^"^ iz) = z'^F{z)Giz). 



Hence, 



(/ * gi; to))^'^^ it) = C:^' [z-F{z)Giz)] (t) = C^' [z^F{z)] (t)£^i [z^G{z)] (t) 
= {f^-/'Ug(-/'\.,to)) it). 
Equivalently, we can write 

[z'^Fiz)G{z)] (t) = C^' [z'^Fiz)] it)C:^' [G{z)] (t) = * g{;to)) it) . 



□ 



5.4 State of the Art 

The results of this chapter are already published in [35] and were presented by the 
author in a seminar with the title Fractional Derivatives and Integrals on Time Scales via 
the Inverse Generalized Laplace Transform on June 9, 2011, at the Faculty of Computer 
Science of Bialystok University of Technology, Bialystok, Poland. 
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Chapter 6 



Conclusions and future work 



In this thesis we define fractional sum and difference operators for the time scales T = Z 
and T = {hZ)a and prove some properties for them. We also established and proved first 
and second order necessary optimality conditions for problems of the calculus of variations 
on both time scales. They show that the solutions of the fractional problems coincide with 
the solutions of the corresponding non-fractional variational problems when the order of 
the discrete derivatives is an integer value. We believe that the present work will potentiate 
further research in the fractional calculus of variations where much remains to be done. 
With respect to the fractional derivatives and integrals on an arbitrary time scale T, we 
introduce a fractional calculus on time scales using the theory of delta dynamic equations 
and the inverse Laplace transform on time scales. The basic notions of fractional order 
integral and fractional order derivative on an arbitrary time scale are given, and some 
properties of the new introduced fractional operators are proved. 

Some possible directions of research are: 

• to extend the results on fractional variational problems for higher-order problems 
of the calculus of variations |27 t H H l52 | l8l ] with fractional discrete derivatives of any 
order; 

• to explore new possibilities of research for fractional continuous variational problems, 
in particular, to get a fractional continuous Legendre necessary optimality condition 
(an interesting open question); 

• to explore several areas of application, for example, in signal processing, where frac- 
tional derivatives of a discrete-time signal are particularly useful to describe noise 
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processes |99| : 

• to find a general formulation leading to the specification of ARMA parameters (as 
pointed out by one of the referees of |33|); 

• to study optimality conditions for more general variable endpoint variational prob- 
lems [65 | [79 | l80 ] and isoperimetric problems [T2 | [T3 ] : 

• to obtain fractional sufficient optimality conditions of Jacobi type and a version of 
Noether's theorem for discrete-time fractional variational problems; 

• to explore direct methods of optimization for absolute extrema |38 [ [5T |[TT3] : 

• to generalize our first and second order optimality conditions to a fractional La- 
grangian possessing delay terms |1]|29]: 

• to extend our results of discrete fractional variational problems to problems involving 
Caputo difference operators [T^ or even Riemann-Liouvile and Caputo together on 
the same Lagrangian. 

Another line of research, where we just have few things done [36], is related with the 
new discrete fractional diamond sum (cf. Definition 11421 below) in the following topics: 

• to investigate the discrete fractional diamond difference introduced in |36] and gen- 
eralize the new discrete diamond fractional operator to an arbitrary time scale T; 

• to investigate the usefulness of modeling with fractional equations |1U1] and study 
corresponding fractional diamond variational principles; 

• to extend the results presented in Chapters [3] and H] to an arbitrary time scale; 

• to study and develop results on other subjects of mathematics using the theory of 
discrete fractional calculus. 

All the results in Chapters [3], H] and [5] were obtained for delta-calculus but could easily 
be rewritten using the operator nabla. However, in our opinion, doesn't make sense to have 
results for A and for V separately. Motivated by the diamond-alpha dynamic derivative on 
time scales |16 [ l77 t l5B] and the fractional derivative of [75 | [55 ] . we introduce a new operator 
that unify and extend the operators A and V by using a convex linear combination of 
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them. To exemplify what we propose, we introduce a more general fractional sum operator 
(for T = Z), making use of the symbol, (cf. Definition I142p . With this new operator 
it's not our intention to change what was done in A and V CcLSGS clS particular cases, 
results on delta and nabla discrete fractional calculus are obtained from the new operator 
- but just to unify the two cases in one. 

Looking to the literature of discrete fractional difference operators, two approaches 
are found (see, e.g., |22|[33]): one using the A point of view (sometimes called the for- 
ward fractional difference approach), another using the V perspective (sometimes called 
the backward fractional difference approach). When 7 = 1 our operator is reduced to 
the A one; when 7 = the operator coincides with the corresponding V fractional sum. 

Analogously to Definition HI one considers the discrete nabla fractional sum operator: 
Definition 139. [22] The discrete nabla fractional sum operator is defined by 



where (3 > 0. Here / is defined for s = a mod (1) and V f is defined for t = a mod (1). 

Remark 140. Let = {a, a + 1, a + 2, . . .}. The operator V^^^ maps functions defined on 
Na to functions defined on Nq. The fact that / and V~^/ have the same domain, while / 
and A^"/ do not, explains why some authors prefer the nabla approach. 

The next result gives a relation between the delta fractional sum and the nabla fractional 
sum operators. 

Lemma 141. J2^ Let Q < m — 1 < v < m, where m denotes an integer. Let a be 
a positive integer, and y(t) be defined on t & Na = {a, a + l,a + 2, . . .}. The following 
statement holds: {A^^y) {t + u) = iV^^y) (t), t G N^. 

With the help of this lemma we are able to rewrite our results of Chapter [3] to the 
V operator. Our intention is not to have results for A and V separately but, instead, to 
extend and unify them through the following operator: 

Definition 142. The diamond-7 fractional operator of order (a, /3) is given, when applied 
to a function / at point t, by 




1 



S = i 



a 



{.Oa'-'-'f) (t) = 7 (Ar/) (t + «) + (1 - 7) i^a'f) (t), 
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where a > 0, /3 > 0, and 7 G [0, 1]. Here, both / and 7<)a"' ^/ are defined for t = a 
mod (1). 

Remark 143. Similarly to the nabla fractional operator, our operator 7<)a"'~'^ maps func- 
tions defined on to functions defined on Na, Na = {a, a + 1, a + 2, . . .} for a a given real 
number. 

Remark 144. The new diamond fractional operator of Definition 11421 gives, as particular 
cases, the operator of Definition |2] for 7 = 1, 

{i^~'''-^f){t) = {A~y){t + a), t^a mod(l), 

and the operator of Definition 11391 for 7 = 0, 

{oOa''^-^f){t) = {V-^f){t), t^a mod(l). 

The next theorems give important properties of the new, more general, discrete frac- 
tional operator 7<)a"'~^. 

Theorem 145. Let f and g be real functions defined on N^, Na = {a, a + 1, a + 2, . . .} for 
a a given real number. The following equality holds: 

Proof. The intended equality follows from the definition of diamond-7 fractional sum of 
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order (a, 

UOa'-^-'if + 9)m = 7(Ar (/ + 9))it + «) + (!- 7)(V:^(/ + 9)m 



1 



1-7 

m 



Y.^t-p{s)Y~\f{s)+g{s)) 



1-7 



7 



□ 



Theorem 146. Let f{t) = k on Nq, k a constant. The following equality holds: 



r(a + l)r(t-a + l) ' "r(/3 + i)r(t-a+i) 

Proof. By definition of diamond-7 fractional sum of order (a, we have 



7 



r(a 

7 
7 



s=0 

r(t-a + l + Q;) 
ar(a)r(t-a + l) 
r(t - a + 1 + a) 
r(a + l)r(t-a + l) 



s=0 



I n r(t-a + i + /j) 
^ + ^^/3r(/3)r(t-a + i)^ 
I n r(t-a + i + /3) 
^ + ^^r(^ + i)r(t-a + i)^- 



Corollary 147. Let f[t) = k for a certain constant k. Then, 

T{t - a + 1 + a) 



(Ar/)(i + «) 



r(a + l)r(t-a + l) 



k. 



□ 



(6.1) 
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□ 



Proof. The result follows from Theorem 11461 choosing 7 = 1 and recalling that 

□ 

Remark 148. In the particular case when a = 0, equality (16. ip coincides with the result 
of |86l Section 5]. 

The fractional nabla result analogous to Corollary 11471 is easily obtained: 
Corollary 149. If k is a constant, then 

(V-^k)(t) - _Jl^^±l+^j, 

Proof. The result follows from Theorem 11461 choosing 7 = and recalling that 

(oOr'^^fc)(t) = (V;^A;)(t). 

Theorem 150. Let f be a real valued function and ai, a2, Pi, P2 > 0. Then, 

+ (1-7) (,0-(-+/'^).-('^^+/'^)/) (t). 
Proof. Direct calculations show the intended relation: 

+ (i-7)(v;^n7^r^'-'v))w 

= f (Ar^ (ArV))(t + «i + «2) + 7(1 - 7) (A;- {V-/^f)) (t + a,) 

+ (1 - 7)7 (V;^^ (A;-/)) (t + «2) + (1 - if (V;^V)) (t) 

= 7' (A;("^+"^)/) (t + ai + a^) + 7(1 - 7) (A,-"^ (A^^V)) (t + «i + /^s) 

+ (1 - 7)7 (v:^^ (vrv)) (t) + (1 - 7)' (v,-(^^+^^)/) (t) 

= 7' (A;(-^+°^)/) (t + «i + «2) + 7(1 - 7) (A-("^+^^V) (t + «i + /32) 

+ (1 - 7)7 (v:(^^+"^)/) (t) + (1 - if (v„-(^^+^^V) (t) 

= 7 [7 (A,-(°^+°^V) (t + «i + «2) + (1 - 7) (V,-(^^+"^V) (t)] 
+ (1 - 7) [7 (A;("^+^^V) (t + «i + + (1 - 7) (V;(^^+^^)/) (t)] . 

□ 
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Remark 151. If 7 = 0, then (^❖""^'"'^^(^❖a"''"^' 



/)) {t) 



f) it)- 



Remark 152. If 7 = 1, then (^❖""^•-^^(^❖r''"'''/)) W = (Aa^"'+"'V) (t + «! + aa). 
Remark 153. If ai = 02 = a and /3i = /32 = P, then 

We now prove a general Leibniz formula. 

Theorem 154 (Leibniz formula). Let f and g be real valued functions, < a, /3 < 1. For 
all t such that t = a mod (1), the following equality holds: 

{.Oa'^'-'ifg)) it) = if^ f7) [(^'^) ^^)] ■ + « + fc)] 



fc=0 



where 



(1 - ^) E ( [(^'^) (^)] [(^."^'^'^/) + /3 + ^)] , (6.2) 
r(n + 1) 



r(t;+ l)r(n-i; + 1)' 
Proof. By definition of the diamond fractional sum, 

{.Oa'-'-'ifg)) (t) = 7 (Ar (/^?)) (t + «) + (1 - 7) (v:^(/^)) (t) 



^X^(t + a-a(.))("-i)/(.)^?(s) 



+ 



1-7 

r(/3) 



J2(t-p{s)r-'f{s)g{s). 



By Taylor's expansion of g{s) 



fc=0 



fc=0 



Substituting the Taylor series of g{s) at t, 



7 



r(«: 



J2it + a-a{s)Y-'^fis) 



J2i-^)'it-sY'^ 



.k=0 



i'^'9)it) 
k\ 



1-7 

m 



J2{t-pis)r^'f{s) 



.fc=0 
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Since 

{t + a- - s)^'^ = {t + a- 

{t - pis)f^{t - sY^^ = {t + l3- a{s)Y^+''+^\ 

t 

and (t-sY''^ 

= 0, we have 



s=t-k+l 



^ ' k=0 ' s=a 



OO 7fc \ / :\ t — k 



r(-a + i)r(«) _ r(-/3 + i)r(/3) 



r(-a + A: + l)r(A; + a) r(-/3 + A; + l)r(A; + /3) 
and k\ = r{k + 1), the above expression becomes 



OO 

7 



fc=0 ^ ^ 

+ (1 - ^) E (t^) (V'=^?)(t)(A;(^+'=)/)(t + /3 + fc). 

t._n V / 



□ 



Remark 155. Choosing 7 = in our Leibniz formula (16.21) . we obtain that 
{^a'{f9)m = £ [iV'gm] [(A,"(^+^)/)(t + (3 + k)]. 

k=0 ^ ^ 

Remark 156. Choosing 7 = 1 in our Leibniz formula (16.21) . we obtain that 

(Ar + «) = E [(^'^)(^)] [(Aa-^"+'^V)(t + a + k)]. (6.3) 
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As a particular case of f)6.3|) . let a = 0. Then, recalling Lemma [TUl we obtain the Leibniz 
formulas of |23] . 

The preliminary results presented here, about this new fractional operator, are pub- 
lished in jSn] and were presented by the author at the NSC'IO, 3rd Conference on Nonlinear 
Science and Complexity, Ankara, Turkey, July 28-31, 2010. 
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Appendix A 

Maxima code used in Chapter 3 



The following Maxima code implements Theorem [Sll Examples illustrating the use of 
our procedure extremal below are found in Section 13.41 



kill (all) $ ratprint : f alse$ simpsuin:true$ tlimswitch:true$ 

sigma(t) :=t+l$ 

rho(t) :=t-l$ 

rho2 (t) : =rho (rho (t) ) $ 

Delta(exp,t) :=block( def ine(f 12(t) ,exp) , 
return((f 12(sigma(t))-f 12(t))) )$ 

p(x,y) : = (gamma(x+l) )/(gainma(x+l-y) )$ 

SumL(a,t,nu,exp) :=block( def ine(f 1 (x) ,exp) , 

f l(t)+nu/gamma(nu+l)*sum((p(t+nu-signia(r) ,nu-l) ) *f 1 (r) ,r,(a),((t-l))) 
)$ 

SumR(tl ,b,nul ,expl) :=block( def ine(f2(x) ,expl) , 

f2(tl)+nul/gamma(nul+l)*suin((p(s+nul-sigma(tl) ,nul-l))*f2(s) ,s, (tl+1) ,b) 
)$ 
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DeltaL (a2 , t2 , alpha2 , exp2) : =block ( 

[alphal :ratsimp(alpha2) , a: ratsimp(a2) ,t :ratsimp(t2)] , 
define (f 3(x) , exp2) , define (q(x) ,SumL(a,x, 1 -alphal ,f3(x) ) ) , qO : q(o) , 
o4:float(ev(ratsiinp(Delta(qO,o)) , nouns)) , o41 : subst (o=t , o4) , 
remfunctionCf ) , renifunction(q) , return(o41) )$ 

DeltaR (t3 , b , alphaS , expS) : =block ( 

[alphal : rats imp (alphaS) ,b : ratsimp(b) ,t :ratsimp(t3)] , 
define(f4(x) ,exp3) , define(ql(o) , (SuinR(x,b, 1-alphal ,f 4(x) ) ) ) , 
qlO:ql(z) , 

o5 : float (ev (rat simp ( -Delta ( (SuinR(x,b, 1-alphal ,f4(x) ) ) ,x) ) , nouns ) ) , 
o51 : subst (x=t , o5) , remfunction(f ) , remfunction(ql) , return(o51) )$ 

EL (exp7 , a , b , alpha? , beta?) : =block ( 

[a: ratsimp(a) ,b:ratsimp(b) , alpha: rat simp (alpha?) , beta :ratsimp (beta?)] , 
define (LL(t ,u, v, w) , exp?) , bl : dif f (LL(t ,u, v, w) ,u) , 
sa: subst ( [u=y (sigma(t) ) , v=DeltaL(a,t ,alpha,y (o) ) , 

w=DeltaR(t,b,beta,y(o))] ,bl) , 

b2:diff (LL(t,u,v,w) ,v) , 

sb : subst ( [t=x,u=y (sigma(x) ) , v=DeltaL(a,x,alpha,y (x) ) , 

w=DeltaR(x,b,beta,y(x))] ,b2) , 
sbl : DeltaR (o ,rho (b) , alpha, sb) , sbll : subst (o=t , sbl) , 
b3:diff (LL(t,u,v,w) ,w) , 

sc : subst ( [t=x,u=y (sigma(x) ) , v=DeltaL(a,x,alpha,y (x) ) , 

w=DeltaR(x,b,beta,y(x) )] ,b3) , 
sc2:DeltaL(a,p2,beta,sc) , sc22 : subst (p2=t , sc2) , return(sa+sbll+sc22) 
)$ 

ELt(exp8, a, b, alphas, beta8,t8) := 

rat simp (subst (t=t8,EL(exp8,a,b,alpha8,beta8) ) )$ 

extremal (L , a , b , A9 , B9 , alpha9 , beta9) : =block ( 
[a:ratsimp(a) ,b:ratsimp(b) , alpha :ratsimp(alpha9) , 
beta:ratsimp(beta9) ,A1 :ratsimp(A9) , Bl :ratsimp(B9)] , 
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eqs:niakelist(ratsimp(ELt(L,a,b,alpha,beta,a+i)) ,i,0,ratsiinp((rho2(b)-a))) , 
vars :makelist (y (ratsimp(a+i) ) , i , 1 ,ratsimp( (rho(b) -a) ) ) , 
Xi: [a] ,Xf : [b] ,Yi: [Al] ,Yf : [Bl] , 

X :niakelist (ratsimp(i) , i , 1 ,ratsimp( (rho(b) -a) ) ) , X : append (Xi ,X,Xf) , 
sols:algsys(subst( [y(a)=Al ,y(b)=Bl] ,eqs) ,vars) , 
Y:niakelist (rhs (sols [1] [i] ) , i , 1 ,ratsimp( (rho (b) -a) ) ) , 
Y: append (Yi,Y,Yf) , 

return (makegamma (rats imp (minf act or ial (makef act (sols[l])))) )$ 
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Appendix B 

Maxima code used in Chapter 4 



The following Maxima code implements Theorem [551 and Theorem [HH 

kill (all) $remf unction(all) $remvalue (all) $ 

fpprec: 100$load("draw")$ 

load ( "eval_string" ) $ 

ratprint : f alse$ 

siinpsum:true$ 

simpproduct : true$ 

tlimswitch : true$ 

sigina(t,h) :=t+h$ 

rho(t,h) :=t-h$ 

rho2 (t , h) : =rho (rho (t , h) , h) $ 

Delta(exp,t ,h) :=block( 

def ine(f 12(t) ,exp) , 

return( (f 12(sigma(t ,h) ) -f 12 (t) ) /h) 

)$ 

p(x ,y ,h) : =h-y* ( ( ( (gamma(x/h+l) ) ) / ( (gainma(x/h+l-y) ) ) ) ) $ 
intHZ (f , var , a9 , b9 , h9) : =block ( 
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[a:ratsiiiip(a9) ,b : ratsinip(b9) ,h:ratsinip(h9)] , 

pi : subst ( var=k*h , f ) , 

res : suni(pl*h, k, a/h, rho(b,h)/h), 

return (res) 

)$ 

SumL(a,t,nu,exp9,hO) :=block( 

define (fl(x) ,exp9) , 

li : ratsimp (a/hO) , 

ls:ratsiiiip((t-hO)/hO) , 

11 : float ( (nu/gamma(nu+l) ) ) , 

(h"nu)*f l(t)+ll*sum((p(t+nu*hO-sigma(r*hO,hO) ,nu-l ,hO) ) *f 1 (r*hO) ,r,li,ls)*hO)$ 

SuniR(t,b,nul,expl,hl) :=block( 

define(f2(x) ,expl) , 

11 : ( (nul/gamma(nul+l) ) ) , 

11 : ratsimp (sigma (t , hi) /hi) , 

Is : ratsimp (b/hl) , 

(hl"nul)*f2(t)+ll*sum((p(c*hl+nul*hl-sigma(t ,hl) ,nul 
-l,hl))*f2(c*hl) ,c,li,ls)*hl)$ 

DeltaL (a2 , t2 , alpha2 , exp2 , h) : =block ( 

[alphal : rat simp (alpha2) , a: ratsimp(a2) ,t : ratsimp (t2)] , 
define (f3(x) ,exp2) , 

define (q(x) ,SumL(a,x, 1 -alphal ,f3(x) ,h) ) , 
qO:q(o) , 

o4: (ev((Delta(qO,o,h)) , nouns)) , 

o41 : subst (o=t , o4) , 

remf unction (f 3) , 

remf unction(q) , 

return (o41) 

)$ 

DeltaR (t3 , b , alpha3 , exp3 , h3) : =block ( 

[alphal : rats imp (alpha3) ,h: ratsimp(h3) ,b : ratsimp (b) ,t : ratsimp (t3)] , 
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def ine (f 4(x) ,exp3) , 

define(ql(x) , (SumR(x,b, 1-alphal ,f4(x) ,h))) , 
q01:ql(o) , 

o5 : - (ev( (Delta(q01 , o ,h) ) .nouns) ) , 

o51 : subst (o=t , o5) , 

remf unction(f 4) , 

remf unction (ql) , 

return(o51))$ 

EL (exp7 , a , b , alpha? , beta? , h7) : =block ( 

[a: rat simp (a) ,b : ratsimp(b) , alpha :ratsimp (alpha?) , beta :ratsiinp (bet 

h:ratsimp(h7)] , 

define(LL(t,u,v,w) ,exp7) , 
bl:diff (LL(t,u,v,w) ,u) , 

sa: subst ( [u=y (sigma(t ,h) ) , v=DeltaL(a,t ,alpha,y (t) ,h) , 

w=DeltaR(t,b,beta,y(t) ,h)] ,bl) , 

b2:diff (LL(t,u,v,w) ,v) , 

sb : subst ( [t=x,u=y (sigma(t ,h) ) , v=DeltaL(a,x,alpha,y (x) ,h) , 

w=DeltaR(x,b,beta,y(x) ,h)] ,b2) , 
sbl :DeltaR(o,rho(b,h) , alpha, sb,h) , 
sbll : subst (o=t , sbl) , 
b3:diff (LL(t,u,v,w) ,w) , 

sc : subst ( [u=y (sigma(t ,h) ) , v=DeltaL(a,x,alpha,y (x) ,h) , 

w=DeltaR(x,b,beta,y(x) ,h)] ,b3) , 
sc2:DeltaL(a,p2,beta,sc,h) , 
sc22 : subst (p2=t , sc2) , 
return (sa+sbll+sc22) 
)$ 

ELt (exp8 , a , b , alphaS , betaS , t8 , h) : = 

rat simp (subst (t=t8 , EL (exp8 , a , b , alphaS , betaS , h) ) ) $ 

extremal2 (L , a , b , A9 , B9 , alpha9 , beta9 , h9) : =block ( 

[a:ratsimp(a) ,b:ratsimp(b) ,h:ratsimp(h9) , alpha :r at s imp (alpha9) , 

beta:ratsimp(beta9) ,A1 :ratsimp(A9) , 
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Bl :ratsiiiip(B9)] , 
eqs rmakelist ( 

(ELt (L , a , b , alpha , beta , a+i*h , h) ) , i , , rat simp ( (rho2 (b , h) -a) /h) 
), 

vars :makelist (y (rat simp (a+i*h) ) , i , 1 ,ratsimp( (rho(b,h) -a) /h) ) , 
Xi: [a] ,Xf : [b] ,Yi: [Al] ,Yf : [Bl] , 

X :makelist (rat simp (i*h) , i , 1 ,ratsimp( (rho(b,h) -a)/h) ) , 

X: append (Xi,X,Xf) , 

eqsOl : subst ( [y(a)=Al ,y(b)=Bl] ,eqs) , 

eqsl : ( (eqsOl) ) , 

return(eqsl) 

)$ 

extremal (L , a , b , A9 , B9 , alpha9 , beta9 , h9) : =block ( 

[a: ratsimp(a) ,b : ratsimp(b) ,h:ratsimp(h9) , alpha :r at s imp (alpha9) , 
beta:ratsimp(beta9) ,A1 :ratsimp(A9) , 
Bl :ratsimp(B9)] , 

eqs rmakelist ( (ELt (L, a, b, alpha, beta, a+i*h,h) ) , i ,0 ,ratsimp( (rho2(b,h) -a) /h) ) , 
vars :makelist (y (rat simp (a+i*h) ) , i , 1 ,ratsimp( (rho(b,h) -a) /h) ) , 
Xi: [a] ,Xf : [b] ,Yi: [Al] ,Yf : [Bl] , 

X :makelist (ratsimp(i*h) , i , 1 ,ratsimp( (rho(b,h) -a)/h) ) , 

X: append (Xi,X,Xf) , 

eqsOl : subst ( [y(a)=Al ,y(b)=Bl] ,eqs) , 

eqsl : ( (eqsOl) ) , 

sols :algsys (eqsl, vars) , 

if length(sols)=l then 

(Y:makelist (rhs(sols [1] [i] ) , i , 1 ,ratsimp( (rho(b,h) -a) /h) ) , 

Y: append (Yi,Y,Yf) 

), 

disp(float(sols)) , 

return (makegamma (rat simp (minf actorial (makef act (sols))))) 
)$ 

LegEqH(L, a, b, alpha, beta, h) : =block( 
mu: 1- alpha. 
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nu: 1-beta, 
duurdiff (L,u,2) , 
duv : dif f (L,u, 1 , v, 1) , 
duw: dif f (L,u, 1 , w, 1) , 
dvvrdiff (L,v,2) , 
dvw: dif f (L, v, 1 , w, 1) , 
dww:diff (L,w,2) , 

parti : (h"2) *subst ( [u=y (sigmaCt ,h) ) , v=DeltaL(a, t , alpha, y(x) ,h) , 

w=DeltaR(t,b,beta,y(x) ,h)] ,duu) , 
part2:2*h" (mu+1) *subst ( [u=y (sigma(t ,h) ) ,v=DeltaL(a,t, alpha, y(x) ,h) , 

w=DeltaR(t ,b,beta,y(x) ,h)] ,duv) , 
partS: (h^mu) *subst ( [u=y(sigma(t ,h) ) ,v=DeltaL(a,t, alpha, y(x) ,h) , 

w=DeltaR(t ,b,beta,y(x) ,h)] ,dvv) , 
part4 : subst (t=sigma(t ,h) ,subst( [u=y (sigma(t ,h) ) ,v=DeltaL(a,t, alpha, y(x) ,h) , 

w=DeltaR (t , b , beta , y (x) , h) ] , dvv) ) * (mu*h"mu-h"mu) "2 , 
parts : intHZ(h"3*subst (t=s , subst ( [u=y (sigmaCt ,h) ) , v=DeltaL(a,t ,alpha,y (x) ,h) , 
w=DeltaR(t ,b,beta,y (x) ,h)] ,dvv) ) * ( (1 /gamma (mu+1) ) *mu* (mu-1) 

*p(s+mu*h-sigma(sigma(t ,h) ,h) ,mu-2 ,h) ) "2 , s , sigma(sigma(t ,h) ,h) ,b,h) , 
part6 : 2*h" (nu+1) *subst ( [u=y (sigma(t ,h) ) ,v=DeltaL(a,t, alpha, y(x) ,h) , 

w=DeltaR(t,b,beta,y(x) ,h)] ,duw)*(nu-l) , 
part7:2*h" (mu+nu) *subst ( [u=y (sigma(t ,h) ) ,v=DeltaL(a,t, alpha, y(x) ,h) , 

w=DeltaR(t,b,beta,y(x) ,h)] ,dvw)*(nu-l) , 
parts : 2*h" (mu+nu) *subst (t=sigma(t ,h) , subst ( [u=y (sigma(t ,h) ) , 

v=DeltaL(a,t , alpha,y (x) ,h) ,w=DeltaR(t,b,beta,y(x) ,h)] , dvw))* (mu-1) , 
part9:h"(2*nu)*subst( [u=y (sigma(t ,h) ) , 

v=DeltaL(a,t , alpha,y (x) ,h) ,w=DeltaR(t,b,beta,y(x) ,h)] ,dww)*(l-nu)"2, 
partl0:h"(2*nu)*subst(t=sigma(t,h) ,subst( [u=y (sigma(t ,h) ) , 

v=DeltaL(a,t , alpha,y (x) ,h) ,w=DeltaR(t,b,beta,y(x) ,h)] ,dww)) , 
part 11 : intHZ(h"3* subst (t=s , subst ( [u=y(sigma(t ,h) ) , v=DeltaL(a, t , alpha, y(x) ,h) , 
w=DeltaR (t , b , beta , y (x) , h) ] , dww) ) * ( ( l/gamma(nu+l ) ) *nu* ( 1 -nu) 
*p(t+nu*h-sigma(s ,h) ,nu-2,h))"2,s,a,t,h) , 
return (part l+part2+part3+part4+part5+part6+part7+part8+part9+partl0+part 11) 
)$ 

LegEqMinH (L , a , b , A , B , alpha , beta , h) : =block ( 
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vec: [] , 

vec : append (vec , [subst( [t=a,y(a)=A] ,ratsinip(LegEqH(L,a,b,alpha,beta,h)))] ) , 
vectTemp: [a] , 

for i:l thru length(X)-3 do( 

tl : ratsiiiip(a+i*h) , 

vectTemp : append (vectTemp, [tl] ) , 

res : subst ( [t=ratsimp(tl)] ,ratsimp(LegEqH(L,a,b,alpha,beta,h) ) ) , 

vec : append (vec , [res] ) 

), 

M: matrix (vectTemp, vec) , 

return (M) 

)$ 

LeqEqMinAllH(L,a,b, A,B,alpha,beta,sols,h) :=block( 
a:ratsimp(a) ,b:ratsimp(b) , 
Validos: [] , 

for i:l thru length(sols) do( 

M : LegEqMinH (L , a , b , A , B , alpha , bet a , h) , 

M:subst([y(a)=A,y(b)=B] ,M) , 

M : subst (makelist (sols [i] [k] , k, 1, length(sols [i] ) ) ,M) , 

Sinai : [] , 

c:0, 

for 1:1 thru matrix_size (M) [2] do ( 
Sinai: append (Sinai, [is(M[2] [1]>=0)]) , 
if is(M[2] [l]>=0)=true then(c:c+l) 
), 

if c=matrix_size (M) [2] then ( 
Validos : append (Validos , [i] ) 
) 

), 

disp (Validos) 

) 

$ 

TESTE(L, a,b, A,B , alpha, bet a, h, sols) : =block( 
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a:ratsiinp(a) ,b:ratsinip(b) , 
indices : [] , 
valor: [] , 

for i:l thru length(sols) do( 
fOO:subst( 

[u=y (sigmaCt ,h) ) ,v=DeltaL(a,t, alpha, y(t) ,h) ,w=DeltaR(t,b,beta,y(t) ,h)] ,L 
), 

rOO:ratsimp(intHZ(fOO,t,a,b,h)) , 
r01:subst(sols[i] ,rOO) , 
r02:subst( [y(a)=A,y(b)=B] ,r01) , 
indices : append(indices , [i] ) , 
valor : append (valor , [r02] ) 
), 

Mat :matrix(indices, valor) , 

return (Mat) 

)$ 
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